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Magnetic  insulation  is  crucial  in  the  operation  of 
large  pulsed  power  systems.  This  phenomena  is  explored  in 
the  following  work.  Particular  attention  is  paid  to 
describing  magnetic  insulation  in  realistic  pulsed  power 
systems.  A  theoretical  model  is  developed  that  allows  one 
to  produce  self  consistent  magnetically  insulated  laminar 
flows  in  perturbed  cylindrical  systems  given  only  the 
electron  density  profile.  The  theory  is  checked  and 
justified  by  detailed  comparisons  with  results  from  a 
2-dimensional  electromagnetic  code  —  MASK. 

The  procedure  followed  in  the  theoretical  development 
is  to  use  the  relativistic  Vlasov  equation.  Ampere’s  law  and 
Gauss’  law,  to  obtain  a  relation  between  the  density  profile 
and  the  velocity  profile  for  insulated  flows.  Given  the 
density  profile  and  the  corresponding  derived  velocity 
profile,  a  self  consistent  flow  solution  is  obtained  by 
means  of  Maxwell’s  equations.  It  is  checked  by  taking  a 
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special  case  —  corresponding  to  no  perturbations  —  wliich 
results  in  the  well  known  Brillouin  flow  theory. 

In  addition,  a  detailed  numerical  study  is  undertaken 
in  which  the  magnetic  insulation  threshold  for  various 
pulsed  power  systems  is  determined.  Particular  emphasis  is 
placed  on  determining  the  magnetic  insulation  threshold  of  a 
pulsed  power  system  employing  a  plasma  erosion  opening 
s w i t  ch . 

The  procedure  employed  in  the  computational  study  is  to 
vary  critical  aspects  of  the  pulsed  power  system  and  then 
note  whether  magnetic  insulation  breaks  down.  The  point  at 
which  mangetic  insulation  breaks  down  —  as  a  function  of 
geometry,  load  impedance,  and  applied  voltage  —  is  the 
magnetic  insulation  threshold  for  the  system.  The 
computational  technique  used  is  a  2-dimensional 
electromagnetic  particle  in  cell  code.  The  results  of  the 
study  are  consistent  with  basic  magnetic  insulation  theories 
in  the  limited  regimes  in  which  they  apply. 

Results  from  the  work  include  magnetic  insulation 
threshold  graphs  for  a  variety  of  systems.  These  graphs 
give  clear  new  insights  into  the  relationships  between  the 
applied  voltages,  perturbation  and  load  impedances,  and 
magnetic  insulation.  In  addition,  the  new  magnetic 
insulation  theory  allows  one  to  theoretically  predict  flow 
solutions  downstream  of  a  perturbation.  Additional  research 
opportunities  are  suggested  in  the  conclusions. 
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CHAPTER  1 
INTRODUCTION 

1.1  Concepts  of  Magnetic  Insulation 

Magnetic  insulation  is  crucial  in  the  operation  of 
large  pulsed  power  systems.  This  applies  especially  to 
systems  with  electric  fields  greater  than  approximately  20 
to  50  MV/m  [1];  because,  explosive  electron  emission  occurs 
from  surfaces  subjected  to  such  fields.  These  electrons 
will  short  out  the  system  at  a  rate  consistent  with  the 
speed  of  light  without  magnetic  insulation.  Therefore, 
insuring  magnetic  insulation  in  large  pulsed  power  systems 
is  a  necessity  of  those  in  the  field  and  a  topic  of  this 
work . 

The  essence  of  magnetic  insulation  is  to  externally 
apply  or  se  1  f -cons i s tent ly  generate  magnetic  fields  that  act 
perpendicular  to  the  electric  fields  found  in  the  system  in 
order  to  confine  the  electrons.  An  example  of  magnetic 
insulation  for  a  se 1 f -cons i s ten t 1 y  generated  magnetic  field 
in  a  coaxial  geometry  is  shown  in  figure  1.  The  current 
flowing  in  the  system  (1,1,1)  generates  a  -B-  field. 

21  C  C  V 

The  electric  field  is  in  the  minus  r  direction.  As  a 
result,  the  electrons  explosively  emitted  from  the  cathode 
are  pulled  by  the  imposed  electric  field  towards  the  anode. 
However,  the  magnetic  field  redirects  the  electrons  back 
towards  the  cathode  by  means  of  vx^  forces.  Therefore,  the 
electrons  ?x?  drift  in  the  +z  direction  while  being 
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Figure  1.  Electron  Orbits  in  Magnetically  Insulated  Flow. 

The  primary  electron  orbits  are  shown  along  with  the 
directions  of  the  anode  (I  )  ,  cathode  ^tid  space 
charge  currents  in  a  cylindrical  geometry. 
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magnetically  insulated  away  from  the  anode. 

In  figure  1,  various  electron  orbit  types  are  shown 
[2].  There  are  common,  laminar,  curtate  and  prolate  orbits 
represented.  If  the  electrons  are  emitted  from  the  cathode 
into  a  region  in  which  the  fields  are  uniform  in  z  and  the 
electrons  initially  have  zero  velocity,  common  orbits 
result.  However,  if  the  electrons  are  born  in  regions  of 
non-uniform  fields,  prolate,  curtate  and  laminar  orbits  can 
be  produced.  In  many  pulsed  power  systems,  prolate,  curtate 
and  laminar  orbits  dominate  because  of  inherent  field 
non-uniformities  due  to  unavoidable  perturbations  in  the 
magnetically  insulated  flows. 

So  there  is  no  confusion,  it  is  important  to  point  out 
that  there  are  three  different  types  of  magnetically 
insulated  (cutoff)  flows  referred  to  in  the  literature  [3]: 
1)  se 1 f- 1 imi ted ,  2)  load-limited,  and  3) 

cons  tan t-f lux- 1 imi ted  flows.  An  example  of  self-limited 
magnetic  insulation  is  given  in  figure  2a.  This  figure 
shows  a  power  pulse  propagating  down  a  long  coaxial 
transmission  line.  The  current  passing  from  the  anode  to 
the  cathode  at  the  pulse  front  {R2)  allows  for  the 
insulation  of  the  line  at  Rl. 

An  example  of  load-limited  magnetic  insulation  is  shown 
in  figure  2b.  The  current  passing  through  the  system  is 
limited  to  that  which  passes  through  the  load  which  is 
responsible  for  the  insulation  of  region  Rl.  It  is 
interesting  to  note  that  the  current  passing  through  the 


4 


1 


a . 


INSULATOR 


CATHODE 

I - ILL - 1 - ^ 

Self-limited  magnetically  insulated  flow.  The  current  flow 
from  the  anode  to  the  cathode  at  the  pulse  front  (K2)  allows 
for  insulation  in  Rl. 


INSULATOR 


\  ANODE 

1 

Rl 

J 

CATHODE 

b.  Load-limited  magnetically  insulated  flow.  The  current  flow 
across  R2  (load  or  diode  gap)  results  in  the  insulation  of 
Rl . 


c . 


INSULATOR 


Constant-flux-limited  magnetically  insulated  flow.  A  current 
from  outside  the  system  (I  )  is  forced  through  the  cathode 
and  produces  a  magnetic  field  of  constant  flux  that  insulates 
the  flow  in  Rl. 


Figure  2.  Types  of  Magnetically  Insulated  Flow’s. 
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load  (R2)  of  figure  2b  is  always  greater  than  or  equal  to 
the  current  lost  across  the  pulse  head  (R2)  as  shown  in 
figure  2a  (for  insulated  flows).  This  means  that  the 
electron  flows  in  load-limited  lines  are  at  least  as  tightly 
trapped  as  electron  flows  in  self-limited  flow  lines.  It  is 
also  interesting  to  note  that  if  the  load  impedance  of 
figure  2b  is  increased  to  the  point  that  current  begins  to 
be  excluded  from  the  diode  region  {load,R2),  the  insulation 
in  R1  returns  to  the  self-limited  type. 

Finally,  one  sees  an  example  of  a  constant-flux 
magnetically  insulated  flow  in  figure  2c.  In  this  case,  a 
current  is  forced  through  the  cathode  from  an  outside  power 
source.  This  current  (I^)  produces  a  magnetic  field  that 
remains  constant  and  is  in  the  minus  0  direction.  If  a 
sufficiently  large  voltage  is  impressed  on  the  line  to  cause 
emission  of  electrons,  the  imposed  magnetic  field  will  cause 
an  drift  and  magnetic  insulation.  However,  in  such 

systems,  the  currents  due  to  the  applied  voltage  are  usually 
neg 1 i g i b i e--un 1 i ke  the  previous  two  cases.  The  power  levels 
are  low  and  the  systems  are  almost  always  non- re  1  a t i v i s t i c . 
This  is  why  an  insulating  field  must  be  supplied  by  an 
outside  source. 

In  this  work,  the  large  pulsed  power  systems  of 
interest  [4]  have  relativistic  electron  flows  and 
significant  se 1 f -genera ted  magnetic  fields.  Also,  most  of 
these  systems  attempt  to  insure  that  all  generated  currents 
pass  through  the  load.  For  this  reason,  the  following 
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research  dwells  primarily  on  load-limited  magnetic 
insulation.  To  give  a  little  more  perspective  and  to 
emphasize  the  importance  of  this  work  a  historical  summary, 
outline  of  research,  the  basic  physics  and  the  computational 
tools  involved  are  presented. 

1.2  Historical  Summary  and  Background 

The  theory  of  magnetic  insulation  begins  with  Hull’s 
first  cartesian,  non- r e 1 a t i v i s t i c  single  particle  treatment 
of  magnetron  cutoff  [5].  Relativistic,  self-consistent, 
cartesian  treatments  begin  around  1973  when  Lovelace  and  Ott 
found  a  condition  governing  cons  tan t - f 1 ux  magnetic 
insulation  of  a  diode  [6],  Also,  at  about  this  same  time, 
Ron,  Mondelli,  and  Ros toker ( RMR )  developed  the 
"quas i - laminar  flow  theory"  employing  common  orbits  (see 
figure  1)  in  cartesian  coordinates  [7].  In  the  work  of  RMR, 
it  is  assumed  that  the  total  energy  and  canonical  momentum 
across  the  electron  sheath  are  equal  to  the  total  energy 
(W=0)  and  momentum  (P^=0)  at  the  cathode.  With  this 
assumption,  they  find  the  electric  and  magnetic  fields 
across  the  sheath,  which  are  expressed  nicely  in  terms  of 
elliptic  functions.  A  laminar  flow  theory  for  magnetically 
insulated  transmission  lines  (MITL)  was  presented  by  Creedon 
in  1975  for  cartesian,  cylindrical  and  conical  geometries 
[8].  This  theory  also  assumes  the  total  energy  and 
canonical  momentum  across  the  flow  are  equal  to  zero.  It  is 
commonly  referred  to  as  the  Brillouin  flow  theory.  Wang,  in 
1977,  using  a  simple  transformation  showed  how  the  equations 


governing  Brillouin  flows  could  be  transformed  into 
Laplace’s  equation  [9].  This  allows  one  to  obtain  Brillouin 
type  flow  solutions  with  arbitrary  cathode/anode  cross 
sections.  In  1979,  Mendel  proposed  a  theory  in  cartesian 
coordinates  that  allows  for  arbitrary  canonical  momentum  and 
total  energy  profiles  across  the  electron  sheath  [10]  for 
flows  with  a  variety  of  electron  orbits. 

Many  other  papers  have  been  written  in  the  more 
ma thma t i ca 1 1 y  tractable  cartesian  coordinates  [11-15,24]; 
unfortunately,  a  majority  of  experiments  and  applications 
involving  magnetic  insulation  utilize  cylindrical  electrode 
geometries  [16-23].  For  this  reason,  it  is  important  that 
the  theory  done  in  cartesian  coordinates  be  extended  to 
cylindrical  coordinates  if  possible.  In  addition,  it  is 
desirable  to  carry  out  new  theoretical  developments  in 
cylindrical  coordinates.  Creedon  and  Wang  were  able  to 
extend  the  laminar  Brillouin  flow  theory  to  cylindrical 
coordinates.  Swegle  attempted  to  extend  RMR  ’  s  work  to 
cylindrical  coordinates,  but  found  that  an  equilibrium 
"quas i - 1  am i nar  flow  solution"  did  not  exist  [25].  In  this 
work,  the  theoretical  developments  are  done  primarily  in 
cylindrical  coordinates. 

Not  suprisingly,  the  theoretical  efforts  mentioned 
above  accompany  very  active  experimental  pulsed  power 
programs  around  the  world.  One  pulsed  power  effort  of 
particular  interest  being  actively  pursued  by  the  US  [4], 
USSR  [26],  etc.  is  light  ion  beam  inertial  confinement 
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fusion{ICF).  This  program  requires  ion  beam  pulses  with 

1  -i  2 

energy  densities  of  10  Watts/cm  with  rise  times  of  15 

nanoseconds  and  durations  greater  than  or  equal  to  20 

nanoseconds.  Power  levels  of  this  magnitude  incident  on  a  1 
3 

cm  deu te r i um- t r i t i urn  pellet  can  compress  and  heat  the 

g 

target  to  10  degrees  Kelvin  and  thermonuclear  ignition. 
However,  the  magnetically  insulated  transmission  lines(MITL) 
that  transport  this  power  to  a  focusing  region  typically 
have  electric  fields  well  above  the  20  to  50  MV/m  required 
for  explosive  emission  of  electrons.  If  it  were  not  for  the 
exploitation  of  magnetic  insulation,  explosive  electron 
emission  would  cause  the  entire  system  to  short  out. 

In  addition,  to  produce  ion  beams  with  the  required  15 
nanosecond  rise  times  and  durations  greater  than  20 
nanoseconds  at  the  necessary  power  levels  requires 
sophisticated  switching  techniques.  One  technique  that  has 
been  used  successfully  and  is  undergoing  extensive  research 
at  Sandia  National  Laboratories,  the  Naval  Research 
Laboratories,  and  other  laboratories  and  universities  around 
the  world  is  the  plasma  erosion  opening  switch  (PEGS). 

A  PEGS  can  conduct  mega-Ampere  levels  [27]  of  current 
for  tens  of  nanoseconds  before  switching  the  current  to  a 
load  in  less  than  ten  nanoseconds.  A  PEGS  operates  as 
follows.  First  the  PEGS  plasma  is  injected  through  an 
exterior  anode  onto  an  interior  cathode.  Since  a  plasma 
acts  as  a  conductor  any  current  forced  down  the  anode  will 
short  to  the  cathode  and  back  to  the  power  source  --  forming 
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an  inductive  storage  loop.  When  a  large  pulsed  power  source 
is  connected  to  such  a  system,  current  flows  down  the  anode 
to  the  PEGS  where  it  shorts  to  the  cathode.  However,  when 
the  current  flowing  through  the  plasma  becomes  so  large  that 
the  free  charges  in  the  plasma  cannot  carry  the  current  the 
plasma  begins  to  erode  away  from  the  cathode.  As  the  plasma 
erodes  away  from  the  cathode,  the  impedance  between  the 
plasma  bottom  and  cathode  becomes  greater  than  that  across  a 
downstream  load.  As  a  result,  the  generated  current  is 
switched  to  the  load  (see  figure  15  and  Chapter  4  for  more 
information  on  the  various  phases  of  a  PEGS  opening).  Now. 
a  PEGS ’ s  successful  operation  depends  critically  on  magnetic 
insulation.  In  fact,  a  PEGS  would  not  open  or  if  it  did 
open,  it  would  not  remain  open  long  if  magnetic  insulation 
broke  down. 

Unfortunately,  there  are  many  regimes  in  which  magnetic 
insulation  breaks  down  even  for  simple  unperturbed  MITLs. 
These  regimes  are  complex  functions  of  the  applied  voltages, 
line  and  load  impedances,  and  the  anode  and  vacuum  currents. 
In  addition,  a  variable  impedance  structure  like  a  PEGS 
severely  disrupts  the  magnetically  insulated  electron  flow 
and  often  results  in  significant  current  loss.  This  is 
because  an  apreciable  fraction  of  the  current  and  power 
transported  by  the  transmission  lines  in  the  regimes  of 
interest  is  carried  by  the  magnetically  insulated  electrons. 

In  order  to  avoid  regimes  in  which  magnetic  insulation 
breaks  down  and  to  more  clearly  understand  magnetically 
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insulated  electron  flows,  a  detailed  theoretical  and 
numerical  investigation  has  been  carried  out  in  this  work. 
Some  parts  of  Mendel's  theory  [10]  have  been  extended  to 
cylindrical  coordinates.  In  addition,  a  new  laminar  flow 
theory  in  cylindrical  coordinates  is  presented  that  allows 
for  non-zero  canonical  momentum  and  total  energy  profiles 
across  the  flow.  The  new  theory  includes  the  previous  state 
of  the  art  (the  Brillouin  flow  theory)  —  which  predicts 
magnetically  insulated  electron  flows  with  zero  canonical 
momentum  and  total  energy  throughout  the  flow  —  as  a 
special  case.  A  large  number  of  computer  simulations  have 
been  done  that  clearly  illustrate  how  magnetic  insulation  is 
tied  to  the  applied  voltages,  line  and  load  impedances,  and 
structure  impedances  (i.e.  a  PEGS)  tf  the  system. 

1.3  Outline  and  Accomplishments 

The  following  is  a  summary  of  what  this  dissertation 
contributes  to  the  field.  Chapters  1  and  2  illustrate  the 
need  for  and  the  basic  theory  governing  magnetic  insulation 
in  cylindrical  pulsed  power  systems.  Some  of  Mendel's  work 
is  extended  to  cylindrical  coordinates.  In  comparing 
predicticTis  from  the  resulting  equations  to  MASK  (see 
section  2.3),  it  is  found  that  the  magnetically  insulated 
flow  is  laminar  in  nature,  but  has  a  non-B r i 1 1 ou i n  flow 
structure. 

This  motivates  the  derivation,  in  Chapter  3,  of  a 
general  laminar  flow  theory  (in  cylindrical  coordinates)  of 
which  the  Brillouin  flow  theory  is  a  special  case.  The  new 


theory  is  used  to  predict  magnetically  insulated  electron 
flow  properties  for  non-Br i 1 1 ou i n  laminar  flows.  It  is 
justified  by  detailed  comparisons  with  MASK  runs. 

In  chapter  4,  the  electron  flow  downstream  of  an 
opening  PEGS  (or  PEGS  like  structure)  is  analyzed.  It  is 
noted  that  a  PEGS  launches  a  magnetically  insulated  electron 
flow  by  means  of  an  drift.  In  addition,  a  set  of 

equations  is  generated  via  the  laminar  flow  theory  developed 
in  Chapter  3  that  gives  the  fields  and  properties  of  the 
1  aunched  f 1 ow . 

The  criteria  governing  load-limited  magnetic  insulation 
are  given  in  Chapter  5.  These  are  presented  in  a  graphical 
form  for  perturbed  and  unperturbed  cylindrical  pulsed  power 
transmission  lines.  The  threshold  for  the  magnetic 
insulation  of  a  system  with  a  PEGS  is  given.  It  is 
confirmed  by  favorable  comparisons  of  the  theoretical 
predictions  to  an  experiment  done  at  NRL  [31]. 

1.4  Basic  Phy s i cs 

The  electron  flows  in  large  pulsed  power  systems  are 
often  relativistic  and  co I  I i s i on  1 es s .  This  is  in  fact  true 
of  the  electron  flows  downstream  of  the  PEGSs  and  in  the 
MITLs  being  considered  herein.  Therefore,  the  relativistic 
Vlasov  equation  is  used  as  a  starting  point  for  theoretical 
development.  It  is  given  by  equation  (1)  in  MKS  units.  In 
fact,  all  equations  and  physical  relations  given  in  this 
text  are  expressed  in  the  MKS  system  of  units.  The 
relativistic  Vlasov  equation  for  electrons  is; 
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Here  v  is  the  directional  velocity  of  the  electron,  m  = 

- 1 1  -  1  9 

9.1095x10  kilograms  (electron  mass),  e  =  1.0022x10 

coulombs  (electron  charge),  f(x.p.t)  is  the  electron 

distribution  function.  ?  is  the  electric  field  and  ^  is  the 

magnetic  induction. 

The  ^  and  fields  are  self  consistently  determined  via 
Maxwell’s  equations’- 

V .t  =  0  .  {^) 

(5) 


vx^  =  -ep  /dpvf  +  ?  , 

o-'  2  at 

c 

=  -j/d^p  f  . 


-*  a^ 

vx^  =  • 
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The  studies  to  be  undertaken  are  often  in  equilibrium 
or  quas i -equ i 1 ibr i urn  which  allows  one  to  drop  the 
derivatives  with  respect  to  time  in  equations  (1),  (5).  and 

(7).  In  addition,  the  electron  density  is  of  particular 
interest  in  much  of  this  work.  It  is  related  to  the 
electron  distribution  function  f(x,v,t)  by: 

p(x,t)  =  -e/d^p  f(x.v,t)  (8) 

The  physics  embodied  in  the  relativistic  Lagrangian 
will  be  a  basis  for  much  of  the  following.  The  relativistic 
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Lagrangian  is  : 

L  =  -mc^/-r  +  eip(x.t)  -  ev.^(x,t)  ,  (9) 

where  ip(x,t)  is  the  electrostatic  potential  and  A(x,t)  is 
the  vector  potential. 

1.5  Computational  Tools 

The  primary  computational  tools  to  be  used  are  the 
2-dimensional  electromagnetic  particle  in  cell  (PIC)  code 
MASK  [28]  and  the  transmission  line  code  BERTHA  [29].  Since 
MASK  is  the  basis  of  this  works  computational  study,  it  is 
introduced  first.  MASK  is  referred  to  as  an  electromagnetic 
code  because  it  solves  the  fully  time  dependent  Maxwell's 
equations  (equations  (4)  -  (7))  allowing  one  to  study 
electromagnetic  waves  with  it.  This  is  to  be  contrasted  with 
electrostatic  or  magnetostatic  codes  which  do  not  include 
the  time  derivatives  seen  in  equations  (5)  and  (7)  and  which 
therefore  can  not  be  used  to  study  systems  in  which 
electromagnetic  waves  are  important. 

MASK  solves  Maxwell’s  equations  in  two  dimensions  and 
in  cylindrical  coordinates  by  dividing  the  region  to  be 
modeled  (i.e.  a  cathode/anode  gap)  into  a  computational 
mesh.  Maxwell’s  equations  are  finite  differenced  on  this 
mesh  with  the  electric  and  magnetic  fields  defined  at  each 
mesh  point.  The  source  of  these  fields  are  charged 
particles  existing  in  and  injected  into  the  region  of 
interest,  along  with  the  specified  boundary  field 
constraints  . 

Typically,  each  charged  particle  (macro  particle) 
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represents  10  to  10  electrons  or  ions.  These 
macro-particles  move  through  the  cells  of  the  computational 
mesh  in  accordance  with  the  fields  defined  at  the 
corresponding  mesh  points.  The  fields  are  used  to  advance 
the  particle  position  within  a  computational  cell  (where  the 
phrase  particle  in  cell  code  originates)  and  vice  versa  by 
means  of  an  integration  technique  called  the  leap-frog 
me  thod . 

The  two  first-order  differential  equations  to  be 
seperately  solved  when  using  the  leap  frog  method  are  [30]: 

m  =  p  ,  (10) 

and 

dx 


d  t 


=  V 


(11) 


Initially  one  must  be  given  the  particle  positions  and 
velocities  along  with  the  boundary  conditions.  Using  a 
weighting  scheme  the  particle  positions  and  velocites  are 
converted  into  charge  and  current  densitii.s.  These,  taken 
together  with  the  boundary  conditions,  are  used  in 
conjunction  with  the  finite  differenced  Maxwell’s  equations 
to  update  the  field  quantities.  After  the  fields  are 
obtained  they  are  used  to  define  the  force  on  a  given 
particle  (the  Lorentz  force): 

?(x. t)  =  -e(?+vx^)  .  (12) 

This  force  is  a  function  of  the  mesh  cell  (position)  in 
question  and  time.  It  is  used  to  obtain  an  average  velocity 


for  the  particle  over  the  time  step  (see  equation  (10)), 
which  in  turn  is  used  to  update  the  particle  position  (see 


15 


equation  (11)).  Given  the  new  particle  positions  and 
velocities,  charge  densities  and  current  densities  are 
obtained  which  are  then  substituted  into  Maxwell’s  equations 
to  obtain  the  time  advanced  fields.  One  should  note  that 
the  charge  densities  used  in  this  process  are  corrected  for 
particle  weighting  errors  by  solving  an  electrostatic 
version  of  Poisson’s  equation.  In  addition,  the  charge 
density  on  the  structure  boundaries  is  modified  to  insure 
the  stucture  potential  remains  as  specified. 

This  is  because  fourier  transforms  (FFT)  are  being  used 
in  the  Poisson  solves.  Since  this  is  so,  the  structures 
within  the  mesh  must  be  defined  in  terms  of  the  surface 
charges  that  result  in  the  required  potential  on  that 
surface.  This  surface  charge  density  must  change  with  the 
changing  plasma  to  maintain  a  constant  potential  on  the 
stucture  surface. 

The  required  surface  charges  are  obtained  be  means  of 
the  capacitive  matrix  method  (30).  This  method  requires  the 
precalculation  of  a  capacity  matrix  C.  This  matrix  relates 
the  plasma  induced  electrode  charge  to  the  potential  of  the 
electrode.  The  procedure  one  follows  is  to  first  solve 
Poisson's  equation  with  no  charges  on  the  electrode  points. 
One  then  records  the  difference  between  the  desired  value 
and  the  potential  actually  found.  If  one  then  multiplies 
this  error  by  C  the  desired  charge  at  each  mesh  point  along 
the  electrode  is  given. 

In  the  following.  MASK  results  are  presented  for 
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cy 1 i ndr i ca 1 1 y  symetric  systems  run  with  and  without  charged 
structures.  Two  sample  MASK  meshes  are  given  in  figure  3. 
Mesh  3a  corresponds  to  the  MASK  runs  for  MITLs  without  any 
perturbations.  Mesh  3b  corresponds  to  runs  for  MITLs  with  a 
structure.  The  cross  hatched  areas  are  space  charge  limited 
(SCL)  ion  emitting  regions  (when  ions  are  desired).  The 
cathode  is  a  space  charge  limited  emitter  of  electrons.  A 
space  charge  limited  emitting  surface  refers  to  a  surface 
covered  by  a  thin  plasma  sheath.  This  is  typical  of  many  of 
the  metallic  surfaces  encountered  in  the  pulsed  power 
transmission  lines  encountered  in  this  study.  Any  electric 
field  normal  to  such  a  plasma  covered  surface  will  cause 
charges  to  flow  out  of  the  surface  plasma  in  such  a  way  as 
to  zero  out  the  normal  component  of  the  field  at  the 
surface.  The  amount  of  charge  required  to  zero  the  normal 
component  of  the  electric  field  is  easily  determined  using 
Gauss’s  law.  MASK  injects  this  amount  of  charge  self 
consistently  at  the  emitting  surface  every  cycle. 

Now  in  figures  3a  and  3b.  V  is  the  anode  voltage;  r 

O  SI 

is  the  anode  radius;  r  is  the  cathode  radius;  r  is  the 

c  p 

radial  distance  to  the  bottom  edge  of  the  structure:  D  is 

the  distance  from  the  left  edge  of  the  mesh  to  the  left  edge 

of  the  structure;  Az  is  the  width  of  the  perturbation,  d  is 

g 

the  diode  gap  spacing  and  D  is  the  axial  length  of  the 
system.  The  above  parameters  are  varied  in  a  series  of  over 
100  calculations  and  discussed  in  chapter  5.  The 
calculations  are  used  to  determine  the  magnetic  insulation 
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a.  The  mesh  set-up  for  simple  MITL  MASK  simulations. 


ANODE  V=Vo 


D 

b.  The  mesh  set-up  for  PEGS  MASK  simulations. 


Figure  3.  MASK  meshes. 
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threshold  as  a  function  of  applied  voltage,  geometry,  and 
load  impedance.  They  are  also  used  to  determine  flow 
properties  downstream  of  a  PEGS. 

This  is  done  by  allowing  the  structure  depicted  in 
figure  3b  to  emit  ions.  This  ion  emitting  structure 
simulates  a  PEGS  at  some  instant  in  time.  In  other  words, 
as  the  switch  plasma  erodes  away  from  the  cathode  it  can  be 
approximated  by  a  conducting  block  that  emits  ions  (SCL). 

In  this  way.  one  can  determine  the  flow  properties  and 
magnetic  insulation  threshold  for  a  system  with  a  PEGS 
without  the  expense  of  modeling  the  PEGS  plasma  itself. 

As  a  check  on  the  magnetic  insulation  threshold 
determined  by  MASK  for  a  system  with  a  PEGS,  a  transmission 
line  code  called  BERTHA  utilizes  the  obtained  threshold  to 
simulate  a  Gamble  II  experiment  [31].  The  simulation 
results  are  found  to  correspond  closely  to  the  experimental 
results  which  partially  justifies  the  ca 1 cu 1 a t i ona 1  study. 

BERTHA  (following  taken  from  reference  (29))  follows 
forward  and  backward  moving  current  and  voltage  waves  and 
then  superimposes  these  waves  to  find  the  physically 
measured  voltages  and  currents  at  any  point  in  a  circuit. 

The  current  and  voltage  at  a  given  point  and  at  a  given  time 
are  given  by: 


o 


and 


V  = 


where  is  the  line  impedance  and  amd  are  the 


(14) 
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voltages  of  the  forward  going  and  backwards  going  waves 
respectively.  From  this,  it  is  clear  that  the  forward  going 
power  flow  is  givenby-' 

P  =  IV  =(V2  -  v2)/z^  .  (15) 

When  the  voltage  waves  are  incident  on  an  impedance 
discontinuity,  as  between  the  MITL  and  load  the  incident 
wave  is  partially  reflected  and  partially  transmitted  (see 
figure  •4).  The  reflection  and  transmission  coefficients  are 
g i ven  by : 


r 


and 


(16) 


t 


(17) 


where  is  the  load  impedance.  Any  number  of  reactive  or 
non-reactive  circuit  components  may  be  included  in  a  given 
problem.  For  more  details  on  this,  see  reference  [29]. 
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CHAPTER  2 

BASICS  OF  CYLINDRICAL  MAGNETICALLY  INSULATED  TRANSMISSION 

LINE  THEORY 

Because  a  majority  of  experiments  involving  magnetic 
insulation  utilize  coaxial  cylindrical  geometries,  and 
because  of  the  lack  of  theoretical  development  in 
cylindrical  coordinates,  magnetic  insulation  in  cylindrical 
coordinates  is  considered  here.  Parts  of  Mendel’s  cartesian 
theory  [10]  are  extended  to  cylindrical  coordinates  and 
checked  using  a  pressure  balance  analysis.  In  comparisons 
of  the  theoretical  results  to  MASK  calculations  and  in 
examination  of  particle  orbits  it  is  determined  that  the 
electron  orbits  in  many  MITLs  are  laminar  in  nature.  As  a 
result,  the  only  laminar  electron  flow  theory  available  in 
cylindrical  coordinates  —  the  Brillouin  flow  theory  —  is 
presented  (no  one  has  yet  been  able  to  produce  a 
self-consistent  non-laminar  flow  theory  in  cylindrical 
coordinates).  Unfortunately,  the  laminar  electron  flows  it 
predicts  do  not  correspond  to  those  obtained  in  a  variety  of 
MASK  computer  experiments.  This  motivates  the  development 
of  a  general  theory  for  laminar  flows  in  cylindrical 
coordinates  which  is  presented  in  Chapter  3. 

2.1  Extend  Cartesian  Theory  to  Cylindrical  Coordinates 
The  nature  of  the  problem  is  shown  in  figure  1.  A 
voltage  difference  (in  MV  regime)  is  set  across  a  coaxial 
transmission  line  and  after  a  few  nanoseconds  an  equilibrium 
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magnetically  insulated  state  is  established.  The  electrons 
that  are  emitted  from  the  cathode  in  a  space  charge  limited 
fashion  enter  a  variety  of  orbits  —  common,  curtate, 

prolate,  laminar,  etc.  -  depending  on  non-un i f o r m i t i e s  in 

the  f i e Ids . 

Since  the  net  flow  of  current  inside  the  transmission 
line  is  in  the  -z  direction  (figure  1),  and  there  is 
azimuthal  symmetry  the  magnetic  field  between  the  anode  and 
cathode  will  be  an  r  dependent  field  in  the  0  direction:  ^  = 
B(r)0.  In  addition,  neglecting  end  effects,  the  gradient  of 
the  applied  potential  yields  an  r-depondent  electric  field: 

?  =  E(r)r.  Because  an  equilibrium  state  is  being  considered 
all  time  dependent  terms  are  dropped  at  the  start  of  the 
theoretical  development. 

Now,  applying  Gauss’s  law  (from  equation  (6))  in 
conjunction  with  the  relativistic  Vlasov  equation  (equation 
( 1 ) )  i mp lies: 


v.f  =  vj^(r)  =  -p(r)/e^  =  /dp^/dp^  .  (18) 

o 

2  d  d 

where  >  "^(r)  is  the  potential  as  a 

function  of  r  across  the  transmission  line,  p(r)  is  the 
electron  charge  density,  is  the  permittivity  of  free 

space,  e  is  the  electron  charge,  p^  and  p^  are  the  z  and  r 
components  of  momentum  respectively,  and  f(p^.p^)  is  the 
p^,p^  space  electron  density. 

It  is  also  true  that  if  ^  =  B(r)0  then: 
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,d  A  d  A  .2 


(19) 


V - ,Q  ~  ^dz.  dr  >  '  ' 

Since  there  is  no  electron  flow  loss  in  the  radial  direction 

and  the  electron  orbits  exist  uni f orm 1 yy ihr oughou t  the 

electron  sheath,  there  is  no  net  radial  electron  flow  as  a 

function  of  r.  This  is  because  for  every  electron  moving 

towards  the  cathode  there  is  one  moving  towards  the  anode  in 

every  macroscopic  control  volume.  As  a  result  =  0. 


Theref  ore 


B{r)0  =  -(|7^z)0  . 


(20) 


This  means  Ampere's  law  (equation  (5))  becomes 

vx^  =  v^A  z  =  p[f  .  (21) 

r  o 

Rewriting  (equati.  (21))  in  terms  of  the  distribution 
function  f(Pj,.p^)  gives: 

v^A(r)  =  p^e/dp^/dpj,  ^  .  (22) 

Tf 

It  is  convenient  to  write  equations  (18)  and  (22)  in  terms 
of  the  canonical  momentum  and  total  energy.  The  appropriate 
Jacobian  of  the  transformation  is  given  by: 


det(^)  = 


^r 

^r 

dp 

cz 

dw 

(23) 

^z 

^z 

^Pez 

dw 

z-componen  t 

of  canonical 

momentum  and  w  is 

the  total  energy  of  a  flow  electron.  The  components  of  the 
Jacobian  can  be  found  using  the  relativistic  Lagrangian 
subject  to  the  constraints  on  the  fields  mentioned 
previously: 
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I  2, ,  ,  .2  .2  2.1/2  ,  .  ,  ,  ^ 

L  =  -me  (1  -  {  r  +  z  )/c  )  +  e<f  {  r  ) -eA^  {  r )  .  (24) 

The  z-component  of  canonical  momentum  and  momentum  are 


respec  t i ve 1 y : 

dL 


p  =  —  =  mTZ  -  eA  (r) 
cz  z'^  ' 

az 


(25) 


and 


p^=  m-TZ  =  p^^  +  eA^(r)  .  (26) 

The  r-component  of  canonical  momentum  and  momentum  are  given 
by  : 


c  r 


5  f 


=  mnrr  =  p 


(27) 


The  total  energy  is  easily  shown  to  be: 

w  =  (7-l)mc^  -  e<p(r)  (28) 

For  convenience  these  relations  arc  written  in  dimensionless 

form  in  terms  of : 

2 


W  =  w/mc‘ 

4)  =  l  +  e<p/mc^ 


a  =  eA  /me 
z 

P  =  p  /me 

X 


(29) 

(30) 

(31) 

(32) 


th 


where  the  subscripted  x  of  equation  (32)  refers  to  x 
component  of  momentum. 

In  this  new  notation,  equations  (25)-(2S)  become: 


P  =  P  -  a 
cz  z 


P^r  =  +  4))^  -  1  -  (P^^  +  a)^)^''^ 

w  =  t4)  . 

This  means  the  Jacobian  is  (from  equation  (23)): 


(33) 

(34) 

(35) 


det(|§i)  =  - 


3x  . 
J 


W+(t) 


((W+4))  -l-{P^^  +  a)  ) 


2,  1/2 


(36) 


Equations  (18)  and  (22)  can  now  simply  be  written 
non-dimensionlessly  in  terms  of  the  total  energy  and 


canonical  momentum: 


v^4)  = 

r  ^ 


P  A(r.P^^.W)(4HW)dP^^dW 


((4HW)2-l-(P^^+a)2)^^2 


V  a  = 
r 


-  rr  A(r.P^^.W)(P^^+a)dP_dW 


C2 


CZ 


((4HW)2-l-(P^^  +  a)^)^''^ 


All  that  remains  is  to  define  A(r.P^^.W): 


A(r.P_.W)  H  i  ^|j(p^^.W)|  H  i  F'{Pcz-'^> 


CZ 


me 


CZ 


(37) 


(38) 


(39) 


This  expression  follows  from  arguments  found  in  Appendix  A 
(which  is  a  derivation  of  equations  (37)  and  (38)  from  a 
slightly  different  point  of  view). 

In  order  to  avoid  rewriting  the  large  integral 
arguments  of  (37)  and  (38)  the  following  definition  is 


given: 

G(4).a)  = 

This  allows  one  to  write  (37)  and  (38)  as 


F(P^3.W)(((}HW)2-l-(P^^  +  a)2)^/2^p^^^^, 


^ 


(41) 


and 


.i  acet-a) 

r  r  5a 


Now 


(42) 


multiplying  equation  (41)  by  and  equation  (42)  by 


and  adding  (with  =  ’)  gives: 


[r(4)-^-a’^)]  '  +  (4)-^-a’^)  =  2G  (43) 
Another  useful  equation  in  the  following  work  is  obtained  by 
multiplying  equation  (43)  by  r  and  regrouping: 

[r2(4)’2-a’2)]’  =  2((rG)--G) 


(44) 
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Since  most  of  the  theory  in  this  dissertation  is  related  to 
these  equations  it  is  crucial  that  they  be  true.  In  the 
next  section,  a  pressure  balance  study  is  perfoimed  that 
justifies  equations  (43)  and  (44).  Further  justification  is 
given  in  section  2.4  where  the  well  known  Brillouin  flow 
results  are  obtained  by  taking  a  limiting  case  of  equation 
(43). 

2.2  Pressure  Balance  Check 

In  this  section,  equation  (■43)  is  used  to  obtain  the 
electric  field  and  magnetic  induction  at  the  sheath  edge. 
This  is  compared  to  the  expressions  obtained  for  the  fields 
using  the  Maxwe 1 1 -S t r e s s  Tensor.  The  fields  are  obtained 
from  the  stress  tensor  by  noting  that  the  electron  cloud 
flows  down  the  axis(z)  in  equilibrium.  Since  the  electron 
sheath  neither  moves  toward  the  anode  nor  the  cathode  the 
sum  of  forces  must  be  zero  (see  figure  5): 

^  ~  d t ^ ^mechan i ca 1  ^  ^field^  ~  ^ 

^  ^vol  d  X  =  /  n.f  ds  .  (^5) 

Now,  ^  ^  -  I  f(E^+B^))  .  (46) 

=  Maxwe 1 1 -S t re s s  Tensor 

where  P  refers  to  the  momentum.  If  one  substitutes  the 
fields  ?  =  E(r)r  and  ^  =  B(r)G  into  equation  (46)  and  then 
substitutes  the  result  into  equation  ('^5)  the  following  is 


f  ound  : 


?1 


O  uj 


X 

CL  5J- 
O  UJ 


Figure  S,  MITL  for  Maxwell-Stress  Tensor  Analysis. 
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*  J'^'"[(B(r)^-E(r)2)e| 


0/2 


-  (B(r)2-E(r)2)e 


d  r 


-0/2-^ 


rd0 


(47) 


where  r  =  -isinG  +jcos0  and  0  =  -icos0  -  jsinG.  After 
integrating  (47),  rearranging  terms  and  assuming  space 


charge  limited  emission  (E{r^)  =  0)  one  obtains: 


r-/  ^2  ,2  r  v2,  ^ 

E(r  )  =  c  (B(r  )  -  — c  B(r  )  )  + 

m^  ''m  r  c'  ' 

m 

(^8) 

^  r  c 

c 

Now,  if  one  integrates  equation  (43)  over  the  entire 
sheath  .equation  (48)  is  again  obtained.  This  indicates 
equations  (43)  and  (44)  are  indeed  true  statements. 

It  is  interesting  to  note  that  the  integral  term  on  the 

right  hand  side  of  equation  (48)  is  composed  of  terms  due  to 

electromagnetic  stresses  in  the  flow  (attributable  to 
non-laminar  orbits  and  the  cylindrical  geometry).  This  term 
is  non-existent  in  cartesian  coordinates.  The  expression  on 
the  left  hand  side  in  equation  (48)  and  the  first  term  on 
the  right  are  pure  pressure  terms. 

2.3  Laminar  Nature  of  Electron  Flows 

In  this  section,  equation  (44)  is  integrated  and  the 
magnitudes  of  the  stress  tensor  terms  ( of  equation  (49)) 
are  found  to  be  relatively  small.  In  addition,  particle 
orbits  for  representative  MASK  MITL  calculations  are 
examined.  The  stress  tensor  check  and  the  simulations 
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indicate  very  laminar  electron  orbits.  This  is  exploited  in 
chapter  3. 

First,  integrating  equation  {^'l)  from  the  cathode  (r^,) 
to  the  sheath  edge  yields  (see  figure  5 )  ■- 

-  r^{4)’^-a'^)  =  [  G  dr  =  P 
i  m  m'^  c^c  c  '  J  r 


m 


m 


(49) 


where  the  subscripts  m  and  c  refer  to  the  sheath  edge  and 

cathode  radii  respectively.  Writing  equation  (49)  in 

dimensional  form  and  noting  that  B  =  *^o^a  and  B  =  ^o^c 

m  -  c  - 


2Trr 


2Tr  r 


m  c 

(by  Ampere's  law  where  I  is  the  cathode  current  and  I  is 


the  anode  current)  implies: 

cp 


t  = 


m 


2Trr 


o  ^j2  _  j2jl/2 

3i  C 


m 


47rcP 


1- 


m 


o  'a  c ' 


1/2  . 
r 


(50) 


where  ?  is  the  electric  field  at  the  sheath  edge  (r  ).  It 

is  also  true  that  (from  Gauss’s  law): 

^  I  -I 
E_  a  c  r  . 
m  =  ^ 


2Tre  r  v, 
o  m  d 


(51) 


where  v^  is  the  average  drift  velocity  of  the  electron  flow. 
From  (50)  and  (51)  one  sees  that: 


Vd  =  c 


rl  -I  1 

1/2  [■ 

a  c 

I  +I 
'•a  c-‘ 

4TrcP 


1  - 


m 


1-1/2 


2  2 
U  (I  -I  ) 
o  a  c ' 


(52) 


It  is  also  obvious  that  the  electric  field  at  the  anode  is 
given  by: 

r  F 

(53) 


V  ^  ^ 

L  =  m  m 

a  - 


In  table  1,  the  electric  field  at  the  anode  (E  )  for 

a  ’ 

various  MASK  simulations  corresponding  to  figure  3a  and 
Appendix  Cl  and  C2  are  given.  In  addition,  the  electric 
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field  occuring  at  the  anode  for  the  corresponding 
theoretical  laminar  flow  (using  MASK  generated  anode  (I  ) 

Si 

and  cathode  (1^,)  currents)  is  presented  for  comparison.  The 
laminar  flow  electric  field  at  the  anode  can  be  obtained 
from  equations  (50)  and  (53)  if  one  realizes  of  (50) 
equals  zero  for  laminar  flows  because  the  radial  component 
of  momentum  equals  zero  for  (P^  =  0)  all  laminar  flow 
electrons.  From  equation  (34)  and  (40)  this  implies 
G{4),a)=0.  and  from  (49)  that  P^  =  0. 

The  percentage  difference  between  the  MASK  and  laminar 
electric  fields  at  the  anode  are  given  in  Table  1.  It  is 
clear,  from  the  table,  that  the  electric  fields  are  close  ( 
within  plus  or  minus  5  percent  except  for  simulations 
M407-M412)  to  what  they  would  be  for  laminar  flows.  In 
figures  6a  and  6b  one  can  examine  the  electron  orbits  for  a 
coujjle  of  these  problems  and  see  their  laminar  nature. 


Lam i nar 


S i mu  1  a  t i on  I 

I 

E 

E  % 

Dif  f 

a 

c 

a 

a 

M407-M412 

-1 .56x10^  - 

1 . 46x10^ 

-6.1x10^ 

-6.6x10^ 

7.5 

M417.M419 

-7.86x10^  - 

3. 44x10"^ 

-8.3x10^ 

-8 . 5x10^ 

2.2 

M418,M420 

-8.46xl0'^ 

-6. 5x10"^ 

-6.8x10^ 

-6.5x10^ 

-4.8 

M414 

-9.85x10^ 

-8.4xl0'^ 

-6 .3x10^ 

-6.2x10^ 

-1.6 

M424 

-7.96x10^ 

-S.Sxlo"^ 

-8.5x10^ 

r- 

o 

X 

CO 

1 .0 

M458 

-1 .77x10^^ 

-1.6x10^ 

-9.4x10® 

-9.3x10® 

-1.4 

Tab  1 e  1  . 

Comparison  of 

anode  e 

lectric  fields  from  MASK 

and 

laminar  theory  (I  and  I  have  the  same  sign  because  they 

SI  c 

are  defined  in  terms  of  a  system  current  loop). 
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For  figure  6a  (simulation  H4.24 ,  see  Appendix  Cl),  it  is 
clear  that  many  of  the  emitted  electrons  are  traveling  .15 
to  .2  meters  before  returning  to  the  cathode.  For  common 
orbits  (see  figure  1).  the  electrons  will  travel 

A 

approximately  2  Larmour  radii  in  the  z  direction  before 
returning  to  the  cathode.  For  this  problem,  the  average 
Larmour  radius  is  approximately  r^  =  .0035  m.  This  means 
the  electrons  are  traveling  approximately  60  Larmour  radii 

A 

in  the  z  direction  before  returning  to  the  cathode.  The 
orbits  are  very  laminar  as  a  result.  Perturbations  and 
gradients  in  the  electric  and  magnetic  fields  are 
responsible  for  this  state  of  affairs. 

For  figure  6b  (simulation  M414),  the  electron  orbits 
are  confined  more  closely  to  the  cathode.  This  increases 
the  likelihood  of  an  electron  returning  to  the  cathode. 
Nevertheless,  the  orbits  shown  in  6b  are  curtate  and 
prolate.  They  are  periodic,  with  a  period  in  z  of 
approximately  .01  m.  The  average  Larmour  radius  is 
approximately  .001  m;  so.  the  electrons  travel  10  Larmour 
radii  in  z  per  gyration.  Their  total  displacement  in  r  per 
gyration  is  approximately  1.5  Larmour  radii.  Therefore,  the 
orbits  are  also  very  laminar  in  nature. 

The  laminar  nature  of  the  electron  flow  observed  with 
MASK  in  cylindrical  coordinates  is  consistent  with  MAGIC  (a 
2-dimensional  electromagnetic  PIC  code)  simulations  done  in 
cartesian  coordinates  at  Sandia.  In  the  Sandia  simulations. 
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it  was  found  that  MITL  flows  in  cartesian  coordinates  are 
laminar  in  nature  but  non-Br i 1 1 ou i n  [32].  By  non-Br i 1 1 ou in 
it  is  meant  that  the  canonical  momentum  and  total  energy  are 
nor.-^ero  across  the  electron  sheath.  The  a.s  sump  t  i  on  of  a 
canonical  momentum  and  total  energy  equal  to  zero  across  the 
sheath  is  fundamental  to  the  Brillouin  flow  theory.  In  this 
work,  the  non-Br i 1 1 ou i n  nature  of  the  flow  is  ascertained  in 
cylindrical  coordinates  by  examining  a  laminar  flow  and 
comparing  it  to  the  Brillouin  solution.  All  the  simulations 
presented  in  Table  1  are  non-Br i 1 1 ou i n  laminar  flows.  This 
is  shown  in  detail  for  cases  M'ISS.  M420  and  M424  by 
comparing  their  magnetic  field  profiles  to  a  Brillouin  flow 
solution  in  section  3.4. 

2.4  Brillouin  Flow  Theory  Derivation 

In  this  section  the  Brillouin  flow  theory  is  derived  as 
a  special  case  of  equations  (40)  and  (43).  This  derivation 
is  done  to  familiarize  the  reader  with  the  state  of  the  art 
in  cylindrical  MITL  flow  salutions  which  will  become  a 
special  case  of  the  more  general  theoretical  developments  to 
be  presented  in  Chapter  3. 

For  this  derivation,  it  is  assumed  that  the  radial 
component  of  velocity  (f)  for  electrons  in  the  flow  field  is 
zero  (as  it  must  be  for  laminar  flows).  It  is  also  assumed 
that  the  total  energy  and  canonical  momentum  across  the  flow 
are  equal  to  zero.  This  is  a  characteristic  of  Brillouin 
flows  and  comes  from  assuming  the  electrons  have  zero  energy 
and  momentum  at  the  cathode,  and  energy  and  momentum 
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conservation.  It  is  also  assumed  that  the  fields  are 
uniform  in  the  z  direction.  The  last  assumption  is 
typically  not  satisfied  in  experiments  because  of  end 
effects,  t ran s i t i ons .  instabilities,  cavity  modes,  and  all 
other  perturbing  effects.  The  non-uniform  fields  are 
responsible  for  non-Br i 1 1 ou i n  flows.  However,  in  the 
following,  one  is  considering  a  non-per turbed  laminar 
electron  fl ow . 

Now.  if  one  sets  the  total  energy  and  canonical 
momentum  equal  to  zero  in  equation  (40)  one  obtains: 

G(4).a)  H  u  =  ((t)^-l-a^)^''^  .  (54) 


and 

u’  =  44' -aa’ 

u 

This  implies  equation  (43)  becomes: 

(4)’2-a’2)  +  [r(4)’2-a-2)]-  =  2u  * 

It  is  important  to  note  that 


(55) 


(56) 


r  =  ^  (4)"-l-a^) 


2^1/2 


cu 


This  implies  u  =  0,  which  forces  the  derivative  of  u  with 
respect  to  r  to  also  be  equal  to  zero:  u’  =  0. 

With  u  equal  to  zero,  equation  (54)  becomes: 


4>^-a^  =  1 


(58) 


2  2 

In  addition,  with  f(r)  =  (4)’  -a'  ),  the  integration  of 

equation  (56)  gives: 


f(r)  =  (const/r)^  .  (59) 

Relation  (58)  leads  one  to  try  a  solution  in  which 
4)=coshS(r),  and  a=sinhS(r).  This  coupled  with  (59)  implies 
S’(r)^  =  (const/r)^  or  that: 
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^  ^  ^  _  £.9.D  ^  ^  ^  which  implies  S{r)  =  const2  Inr  (60) 
To  simplify  matters,  space  charge  limited  emission  is 
assumed  so  that  the  electric  field  at  the  cathode  is  zero. 
Also,  the  electrostatic  potential  and  vector  potential  ai e 
conveniently  set  to  zero  at  the  cathode.  The  magnetic  field 
at  the  cathode  is  given  by  Ampere’s  law  (B^  =  p^I^/2Trr^). 

The  total  current  flowing  in  the  system  is  governed  by  the 
potential  at  the  sheath  edge  and  at  the  anode  in  conjunction 
with  the  system  geometry.  In  a  somewhat  lengthj’  derivation 
[8],  Creedon  relates  these  quantities  for  a  Brillouin  flow. 
He  finds  the  anode  current  to  be(this  relation  is  obtained 
independently  of  Creedon  in  Chapter  3): 


I  = 
a 


2ir  mc^ 

a 

m 

p.  ce 
o 

ln(ra/rc)l 

m  f,  r  .  /  2  ,^l/2^  -r  -a 

rTTT  )+_a_JL_  • 

‘  ‘  _  y  L  9  1  /9  J 


,  2  , . 1/2 


He  also  shows  that  the  potential  at  the  sheath  edge  for  a 
Brillouin  flow  is  related  to  the  anode  and  cathode  currents 
by  (this  is  also  independently  obtained  in  Chapter  3): 


^m  ^a'^'^c  ’  (62) 

where  a  is  the  relativistic  factor  at  the  sheath  edge  and 
m  ^ 

a  is  the  relativistic  factor  at  the  anode.  a  is  related  to 
a 


the  potentials  by  a  = 


1  +  eV 


Where  V  is  the  electrostatic 


potential  in  question.  Equation  (62)  coupled  with  equation 
(61)  places  a  restriction  on  the  allowed  ratios  of  anode 
currents  to  cathode  currents  in  the  Brillouin  flow  theory. 
The  theory  will  not  satisfy  Ampere's  law  if  equations  (61) 
and  (62)  are  not  simultaneously  satisfied.  An  example  of 
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the  Brillouin  flow  theory  not  satisfying  Ampere’s  law  is 
shown  in  Chapter  3. 

From  equations  {29)-{32).  (59)-(62)  and  the  fact  that 

4)=^coshS(r)  and  2  =  sinhS(r)  one  obtains  the  following 

self-consistent  MITL  flow  solution  ■■ 

2 

m  c 

<p(r)  =  — ^ - (  cosh  (  A  In  (  r /r  ^ )  )  -  1  )  ,  (63) 


A(r)  =  — ^  -  s i nh ( A  1 n ( r/r^ ) )  z 


(64) 


me. 

^(r)  =  — ^ -  —  s inh{ A ln{ r/r^ ) )  r 


(65) 


B(r)  =  — - - -  cosh(Aln{r/r^) )  0 


(66) 


v(r)  =  c  tanh( A ln( r/r^ ) )  z 


(67) 


2  .2 

m  c  e  A 

p(r)  =  -  — - 1 cosh{  A  In  {  r /r^  )  ) 


e  r 


(68) 


cosh 


r  =  r  exp( - r- 

m  c  ^ '  A 


m 


■)  . 


(69) 


A  = 


ln{ra/rc) 


,  -1  ^  TT  — r 

cosh  -r  +  o  m 

,  2  ,  ,  1/2  " 
(  T  -  1  ) 

m  ' 


(70) 


where  <p  is  the  electrostatic  potential,  A  is  the  vector 
potential,  ?  is  the  electric  field,  ^  is  the  magnetic 
induction,  v  is  the  particle  velocity,  p  is  the  electron 


charge  density,  and  r  is  the  electron  sheath  edge  radius. 

m 
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These  relations  were  first  presented  by  Wang  and 
DiCapua  in  1980  [33]  based  on  Creedon’s  findings  given  in 
equations  (61)  and  (62).  Till  now.  this  set  of  equations 
was  me  only  self  consistent  cylindrical  MITL  flow  solution. 

Unfortunately,  the  Brillouin  flow  solution  does  not 
correspond  to  numerical  simulations  done  around  the  country 
—  which  show  laminar  but  non-Br i 1 1 ou i n  flows.  The 
discrepancy  is  due  to  the  fact  that  the  more  realistic 
2-dimensional  simulations  allow  for  perturbations  in  the 
system  not  taken  into  account  in  the  Brillouin  flow  theory. 
This  deficiency  is  eliminated  by  theoretical  de ve 1 opemen t s 
presented  in  the  next  chapter. 
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CHAPTER  3 

GENERAL  THEORY  OF  LAMINAR  FLOWS 

Until  now,  the  only  relativistic  self-consistent  MITL 
flow  solution  in  cylindrical  coordinates  was  given  by  the 
Brillouin  flow  theory  which  applies  specifically  to  laminar 
electron  flows  with  zero  canonical  momentum  and  total  energy 
across  the  sheath.  This  corresponds  to  an  ideal  system  with 
no  perturbations.  Unfortunately,  all  experimental  systems 
have  perturbations  and  end  effects. 

In  this  chapter,  a  general  theory  for  laminar  MITL 
flows  is  presented  in  which  perturbed  flows  can  be  modeled. 
The  Brillouin  flow  theory  is  found  to  be  one  of  an  infinite 
number  of  possible  laminar  MITL  flow  solutions  predicted  by 
the  more  general  theory.  From  the  general  theory,  it  is 
discovered  that  the  specification  of  a  density  profile 
determines  unique  values  of  the  total  energy  and  momentum 
across  the  flow.  With  this  in  mind  a  density  profile  is 
generated  based  on  equations  (40)  and  (43).  The  profile  is 
then  used  to  generate  the  complete  Brillouin  flow  solution 
via  the  general  theory  of  laminar  flows.  The  canonical 
momentum  and  total  energy  are  found  to  be  zero  (figure  11b) 
as  expected.  This  comprises  an  independent  derivation  of 
the  Brillouin  flow  theory  and  a  check  of  the  general  theory 
of  laminar  flows. 

The  general  theory  is  then  used  to  simulate  several 
MASK  calculations  for  short  MITL’s  in  which  the  density 
profile  across  the  flow  is  approximately  constant  and 
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therefore  non-Br 1 1 1 ou in .  These  are  considered  because  of 
the  ease  of  generating  a  self-consistent  theoretical  set  of 
equations  for  a  constant  density  profile.  It  should  be 
emphasized  at  this  point  however,  that  any  density  profile 
could  be  substituted  into  the  key  equa t i ons ( ( 75 ) - ( 84 ) )  to 
obtain  a  self-consistent  laminar  MITL  flow  solution. 
should  also  be  emphasized  that  the  specification  of  a 
density  profile  is  equivalent  to  specifying  the  canonical 
momentum  and  total  energy  across  the  .low,  or  specifying  the 
velocity  profile,  or  the  magnetic  field  profile,  etc.  In 
other  words,  given  one  type  of  detailed  information  about 
the  flow  the  general  theory  allows  one  to  obtain  all  other 
pertinent  flow  quantities.  Since  it  is  especially 
convenient  to  do  this  using  the  density  profile,  the 
following  work  is  presented  as  if  one  had  knowledge  of  the 
density  profile  across  the  flow. 

Nevertheless,  if  one  had  detailed  knowledge  of  the 
magnetic  fields  or  current  across  the  sheath  (from 
experiments),  the  following  equations  are  still  useful. 

After  coding  them  on  the  computer,  the  density  profile  could 
be  altered  iteratively  until  the  magnetic  field  or  current 
profile  of  interest  across  the  sheath  is  duplicated.  One 
would  then  have  a  complete  flow  solution  given  a  minimum  of 
inf  o  r ma  t i on . 

3.1  Relation  Between  the  Density  and  Velocity  Profiles 

The  general  theory  is  obtained  via  a  surpris'ng 
relation  between  the  density  and  velocity  profile  for 
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relativistic  MITL  flows.  A  specification  of  the  density 
profile  implies  a  particular  velocity  profile  and  vice 
versa. 

This  interesting  result  is  derived  in  Appendix  B  in 
which  the  following  procedure  is  carried  out"  First, 
equation  (43)  is  integrated  from  the  cathode  to  some 
arbitrary  r  within  the  sheath.  The  integrated  result  is 
rearranged  so  that  the  electric  field  {^(r))  is  set  on  the 
left  hand  side  of  the  equation.  Using  Ampere's  law,  the 
magnetic  field  quantities  on  the  right  hand  side  of  the 
equation  are  written  in  terms  of  the  anode  and  cathode 
currents.  Then,  equation  (44)  is  integrated  from  the 
cathode  to  r  to  obtain  another  expression  for  The 

resulting  equations  are  evaluated  at  the  sheath  edge  (r^) 
and  equated.  After  differentiation  and  algebraic 
manipulation  a  relation  between  the  density  and  velocity  is 
obtained.  This  relation  is  (given  in  Appendix  B): 

f  ( r  )  =  ^  ,  (72) 

1-k  g(r)‘^ 

21 

where  1  =  y  ^  ,  (72) 

a  c 


(74) 
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f(r) 


(75) 


c 

where  f(r)  is  a  normalized  weighted  average  of  the  electron 
charge.  It  equals  zero  at  the  cathode  and  smoothly 
increases  until  reaching  the  value  of  one  at  the  sheath  edge 
(not  to  be  confused  with  the  distribution  function 
in  the  Vlasov  equation). 

Simple  calculus  and  algebraic  manipulation  of  equations 
(7l)-(75)  results  in  a  more  clearly  defined  relation  between 
the  velocity  and  density  profiles  (this  is  done  in  Appendix 
B): 


v(r) 


1 


2  2  2 
(I  -  I‘^)f‘^(r) 
'a  ' 


c 


1/2 


1  + 


(i; 


t1/2 


(76) 


c  J 

where  I  is  the  anode  current.  I  is  the  cathode  current, 
a  c 

p(r)  is  the  density  profile  and  f{r)  is  the  normalized 
weighted  density  profile. 

3.2  Laminar  Flow  Equations 

Given  a  densi^^y  profile  and  the  corresponding  velocity 
profile  (equation  (76)),  the  self-consistent  MITL  flow 
solution  is  easily  obtained  via  Maxwell’s  equations  (subject 


to  the  restrictions  of  section  2.1).  The  resulting 
equations  are  collected  below  in  a  convenient  form.  Given 


the  charge  density  (p(r))  and  the  velocity  (from  equation 
(76))  the  current  density  is  : 

T  =  P(r)v(r)  .  (77) 

The  electric  field  is 


j: 


p(  r )  r  dr 


^(r)  = 


(78) 


e  r 
o 


where  ?(r^)  =  0  ,  because  the  cathode  is  a  space  charge 
limited  emitter  of  electrons.  The  corresponding 
electrostatic  potential  is 


<p(r) 


E(r)  dr 


(79) 


where  <p(r^)  is  chosen  to  be  zero.  The  magnetic  induction  is 
simply  written  as 


g(r)  = 


cr 


j(r)  r  dr  + 


2Tr 


(80) 


The  vector  potential  is 


A(r) 


-j: 


B(  r )  dr  z 


(81) 


where  A(r^)  =  0  is  chosen  to  equal  zero  at  the  cathode. 

Now,  from  equations  (26)  and  (35)  the  canonical 
momentum  and  total  energy  across  the  flow  are  given  by: 
^cz^^^  ^  m-Tv(r)  -  eA(r) 
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m  V  (  r  )  i  /  ^ 

,,  2.  .  .  2.1/2  ' 

(1  -  V  (r)/c  ) 


The  total  energy  is 


(82) 


W{r)  =  -me' 


1  - 


/■I  2,  ,  ,  2.1/2 

(1  -  V  (r)/c  ) 


-  e^( r  ) 


(83) 


Also  for  laminar  flows. 

?(r)  =  -v{r)  X  ^(r)  (84) 

Equations  (76).  (78).  and  (81)  must  satisfy  equation 

(84).  For  this  reason  equation  (84)  is  used  to  check  the 
validity  of  laminar  flow  solutions  produced  by  the  general 
theory.  An  example  of  this  kind  is  shown  in  figure  10b. 

If  the  sheath  edge  radius  equation  (75)  is 

unknown  it  is  easily  found  given  I  .  I  ,  the  applied  voltage 

Q  C 

across  the  gap(V  )  and  equations  (76)-(79).  This  is  done 

for  the  constant  density  flow  derived  in  section  3.3.2. 

From  an  experimental  point  of  view,  it  is  much  more 

desirable  to  have  the  MITL  flow  as  a  function  of  I  .  I  and 

a  c 

V  rather  than  I  .  I  and  r  .  This  is  because  the  applied 
a  a  c  m 

VO  1 tage  across  the  diode  is  usually  known;  the  outer  sheath 
edge  radius  of  the  electron  flow  is  not. 

In  the  foregoing,  it  was  stated  that  given  only  a 
density  profile,  a  self-consistent  MITL  flow  solution  is 
obtained.  This  is  indeed  true;  however,  if  one  desires  the 
particulars  of  a  solution  for  a  given  voltage,  line 
impedance,  and  load  impedance  (or  equivalently  voltage,  line 
impedance,  and  anode  and  cathode  currents)  the  latter  must 
be  supplied.  Equations  (71)-(84)  were  presented  with  this 
in  mi nd . 
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3.3  Two  Derivations 

In  order  to  demonstrate  the  general  theory  of  laminar 
flows,  a  couple  of  interesting  MITL  flows  are  derived. 

First,  the  Brillouin  flow  theory  is  derived  independent  of 
Creedon’s  total  current  relation  (equation  (61)).  The  two 
Independently  derived  theories  are  plotted  against  each 
other  and  are  found  to  match. 

Then,  the  general  theory  is  used  to  generate  a 
self-consistent  MITL  flow  solution  for  flows  with  constant 
density  profiles.  This  is  done  because  in  simulations  of 
short  MIT's,  quas i -cons  tan t  density  profiles  were  observed. 
The  theory  is  then  compared  to  the  simulations  and  found  to 
accurately  predict  quantities  of  interest. 

3.3.1  Brillouin  Flow 

In  this  section  an  independent  derivation  of  the 
Brillouin  flow  theory  is  presented.  Starting  from  equations 
(40)  and  (43),  one  quickly  obtains  equations  (5S)-(60),  and 
the  relations  <})=coshS(r)  and  a  =  sinhS(r)  as  shown  in  section 
2.4.  From  these  and  Maxwell’s  equations  it  is  easy  to  show 
that  the  Brillouin  density  profile  will  have  the  following 
f  orm  : 

2  .2 

me  e  A 

p(r)  =  -  - —2.  cosh(A^ln(r/r^)  )  (85) 

r 

This  is  equivalent  to  equation  (68)  except  that  the  critical 
constant  ”A^"  (equation  (70)),  previously  obtained  from 
Creedon’s  analysis,  is  indeterminate.  Since  Creedon’s 
results  are  not  needed  and  in  fact  are  derivable  from  the 
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general  theory,  the  following  constitutes  an  independent 
derivation  of  the  Brillouin  flow  theory. 

From  equation  (85).  equation  (75)  becomes; 


sinh(A^ln(r/r^)) 

~  sinh(A  ln(r  /r  )) 


n  '  c 
^n  '  m  ■  c 

This  means  the  velocity  (equation  (76))  becomes; 


(86) 


v(r)  = 


|I 


2irmc  £  A  sinh(A  ln(r/r  )) 
on  '  n  c’  ’ 

(l+m,sinh^(A  In—  )) 

1  n  r 

c 


(l+m,sinh^(A  In—  )) 
'  1  n  r  ^  ' 


1/2 


(87) 


2  2 
I  - 

where  m^  =  - - 


I  sinh(A  ln(r  /r  )) 
c  'n  'mc^'' 


(88) 


Now,  m^  must  equal  1.  This  can  be  seen  by  first  noting  that 

the  density,  electric  field,  and  potential  do  not  depend  on 

.  So,  from  the  potential  (equation  (91))  evaluated  at  the 

sheath  edge  and  the  realization  that  i  =1/1  (from 

mac 

equation  (124))  at  the  sheath  edge  it  is  possible  to  obtain 
A^.  One  can  then  solve  equation  (88)  for  A^ .  By  comparing 
the  resulting  equations  it  is  evident  that  m^  equals  one. 
This  implies  equations  (77)  and  (87)  become  (when  r^  <  r  < 

r  )  = 

m  ■' 


el 


v(r)  = 


27rm  c  £  A 
o  on 


tanh(A  ln(r/r  )  z 


(89) 


J(r)  =  - 


A  ll  I 


2tt  r  ‘ 


sinh(A  ln( r/r  )  )  z 


n 


(90) 


From  (78)  thru  (81)  in  the  region  ^  ^  ^ 


after  algebra 
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one  obtains: 


m  c 

«p{r)  =  — - {cosh{A^ln{r/r^)  )-l ) 


(91) 


r(r)  =  - 


2  A 

m  c  A 

- -  ■“  s inh( A^ln ( r/r ^ ) )  r 


(92) 


and 


^  "^7a -  slnh(A^ln(r/r^) )  z 

n 


^  '  “27r -  cosh(A^ln(r/r^)  0 


(93) 


(94) 


In  the  region  r  <  r  <  r  the  following  relations  apply: 
m  “  ~  a 


'p{r)  = 


m  c 
o 


A  ln(r/r  ) 
n  m ' 


2  2 

r  I  -  I  t1/2  I  -  I 
a  c  I  ^  a  c 

+  • 


r 


I 


.  (95) 


?(r)  = 


m  c 

— -  A  sinh(A  ln(r  /r  )) 

er  n  'n  'm  c 


(96) 


U  I  1  I 

^(r)  =  ^  (  1 1^  I  ln(r/r^)  +  -^s  i  nh  ( 1  n  (  r^/ r  ^  )  )  )  .  (97) 


=  -  2^  I'al  ®  ■ 


(98) 


and 


p(r)  =  v(r)  =  j(r)  =  0 


(99) 


where  from  (88)  and  (94)  the  constant  is 


A  = 


In ( r  /r  ) 
m  c 


s  i  nh 


-1 


l2  -  l2 
a  c 


1/2 
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If  one  carefully  compares  the  equations  in  this  section 
to  those  in  section  2.4,  it  can  be  seen  that  the  formulation 
and  source  of  the  constants  is  different  but  the  functional 
forms  are  the  same.  If  one  compares  the  flow  quantities 
predicted  by  these  two  independent  formulations  it  becomes 
apparent  that  they  are  predicting  the  same  flow  quantities 
given  the  same  inputs.  This  is  done  for  the  potentials. 
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fields,  density,  momentum  and  energy  profiles  in  figures  7 
thru  11.  In  the  figures,  BRI  TH  corresponds  to  Creedon’s 
solution  and  BRIN  TH  corresponds  to  this  works  independent 
derivation.  Also  in  figures  7  thru  11  are  the  constant 
density  (RHO  TH)  solution  (derived  in  3.3.2)  results  in 
which  the  anode  and  cathode  currents  along  with  the  total 
charge,  are  the  same  as  the  Brillouin  flow  case.  It  is 
clear  from  the  figures  that  the  two  independent  Brillouin 
flow  solutions  overlay  each  other.  The  contrast  between  the 
constant  density  flow  solution  and  the  Brillouin  flow  is 
apparent  in  these  figures. 

The  Brillouin  and  constant  density  flow  solutions  shown 
in  figures  7  thru  11  are  characterized  by  the  following: 

=  -  7.18  X  10^  Amperes 

4 

-  -  4.116  X  iO  Amperes 

r  =  .05  Meters 
a 

r  =  . 025  Me  ter  s 
c 

The  Brillouin  flow  cases  have  a  potential  difference  of  = 
0 

2  X  10  volts  between  the  anode  and  cathode  (figure  8a). 
However,  the  constant  density  solution  allows  a  lower 
potential  for  the  same  currents  (figure  8a)  in  this  case. 

Eventhough  the  constant  density  and  Brillouin  flow 
solutions  have  differing  potentials  applied  across  their 
respective  anode/cathode  gaps,  the  anode  electric  fields 


49 


This  relation  is  in  accordance  with  the  laminar  (P  =  0) 
form  of  equation  (50).  Since  the  electric  fields  are  the 
same  at  the  anode  and  the  fields  of  the  Brillouin  solution 
are  larger  close  to  the  cathode  (figure  9a),  the  potential 
for  the  constant  density  flow  is  lower  (figure  8a)  than  the 
Brillouin  flow  for  this  case. 

Before  proceeding  to  the  next  section,  it  is  worth 
discussing  figures  lOb,  11a  and  lib.  Figure  10b  gives  ^  for 
a  Brillouin  flow  field  and  then  compares  it  to  -v  x 
These  overlap  up  to  the  sheath  edge  where  v  goes  to  zero. 
This  verifies  that  the  flow  is  laminar.  This  test  was 

performed  and  passed  on  all  the  laminar  flow  solutions. 

2  2 

Figure  11a  gives  t  curves  for  the  constant  density 

2  2 

and  Brillouin  flow  cases.  nr  -P  must  equal  1  .  as  long  as 

cz 

one  is  within  the  electron  sheath  and  is  considering  a 
laminar  flow.  This  relation  is  seen  to  be  satisfied  in 
figure  11a.  Figure  lib  graphs  the  total  energy  vs  the 
canonical  momentum  for  the  constant  density  and  Brillouin 
flows.  The  constant  density  flows  have  a  profile  in  the 
fourth  quadrant.  The  two  Brillouin  flows  have  canonical 
momentums  and  energies  equal  to  zero  across  the  flow.  The 
reason  the  canonical  momentum  and  total  energy  are  not 
exactly  zero  for  the  Brillouin  flow  solution  obtained  in 
this  section  is  becau  e  insufficient  iterations  (three 
iterations)  were  taken  in  evaluating  equations  (101)-(103). 
The  more  iterations  taken  the  closer  the  canonical  momentum 
and  total  energy  are  to  zero. 


1.6 


a.  Elertrostatic  potential.  ( 


b.  Vector  potential. 

Figure  8.  Potential  Profiles  f< 
Flow  Solutions. 
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b.  Total  energy  vs  canonical  momentum. 

Figure  1  1  .  Laminar  Flow  Check  and  Total  Energy  Vs  Canonical 
Momen  i  Profiles. 
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3.3.2  Constant  Density  Flow 

In  various  MASK  calculations  of  short  MITLs,  it  has 
been  noted  that  the  density  as  a  function  of  radius  is 
relatively  constant  out  to  the  sheath  edge,  where  it  drops 
to  zero.  For  this  reason  (and  the  fact  that  the  integrals 
are  simpler  to  evaluate),  a  self-consistent  MITL  flow  for  a 
constant  density  profile  is  given  below.  Substituting  a 
constant  density  profile  into  relations  (7.5)-(S‘l)  and 
evaluating  between  the  limits  r  1  r  <  r  results  in: 


f{r)  = 


r  -  r 


r  -  r 
m  ( 


-  I 


v(r)  = 


,2  2, 

Trp(r  -r^) 


1  +  (-A-  -)f  (r) 
c 


1/2 


1  +  (r) 


1/2 


jCr)  =  -|p|v(r) 


I  p  I  ,2  2  „  2, 

=  (r  -  r  -  2rln 


'P{r)  =  ^  { 

o 


c  c 


2e  r 
o 


A{r)  =  - 


'"o^c  ,,,  ,  ^ 


nd 


(105) 

(106) 

(107) 

(108) 

(109) 

(110) 
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M  I  rr  I^-I^  T  r  r^-r^  i2 
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where 
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Evaluating  these  same  quantities  from  the  sheath  edge  to  the 

anode  (r  <  r  <  r  )  results  in: 

^  m  ~  a.' 
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and 


v(r)  =  p(r)  =  j(r)  =  0 


(119) 


So,  is  obtained  iteratively  given  the  anode  and  cathode 

currents,  along  with  the  anode  and  cathode  radii  ,  and  the 

voltage  difference  across  the  diode.  Once  r  is  known  all 

m 

other  quantities  defined  in  equations  (105)-(121)  are  given. 

In  typical  experiments  r  ,  r  ,  I  ,  and  V  are  known.  I  is 

a  c  a  a  c 

a  little  harder  to  obtain,  but  can  be  found  using  a  magnetic 


probe  next  to  the  cathode[34]. 

It  is  also  interesting  to  note  that  t  =1/1  at  the 

mac 

sheath  edge.  This  can  be  seen  by  evaluating  (106)  at  the 
sheath  edge: 

=  f  ^  if  -  (123) 

a 

Equation  (123)  gives  the  velocity  of  the  electrons  at  the 

sheath  edge,  so  the  relativistic  factor  at  r  is 

m 


^  =  _ I _  .  la 

(1  -  v2(r^)/c2)I/2 


(124) 


Equation  (124)  also  holds  for  the  Brillouin  laminar 

flow  solution  (see  equation  (62)).  What  thic  means  is  that 

even  though  the  electron  sheaths  have  differing  electron 

density  profiles  and  sheath  edge  radii(r  ),  tlic  electron 

m 

velocity  at  the  sheath  edge  is  the  same  given  the  same  ratio 

of  I  /I 
a  c 

Since  t  =  I  /I  at  the  sheath  edge  for  two  very 
mac  ° 

different  density  profiles,  perhaps  t  =1/1  at  the  sheath 

mac 

edge  for  all  laminar  flow  solutions.  This  is  in  fact  true 
and  is  easily  proved  using  the  general  theory  of  laminar 
flows.  First  one  needs  to  note  that  f(r  )  =  1  at  the  sheath 
edge  (see  equation  (75)).  In  examining  equation  (76)  one 


needs  to  realize  that  the  integral  27r/  p(r)rdr  is  simply 

c 

the  total  charge  (q)  in  the  sheath  per  unit  length  in  z. 


From  equation  (50)  with  P  =0  and  Gauss's  law  the  total 

m 

electron  sheath  space  charge  per  unit  length  is  given  by: 
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(l2  _  j2jl/2 

a  c 


(125) 


Substituting  f(r^)  =  1  and  equation  (125)  into  equation  (76) 
results  in: 


(  \  c  ,.2  t2>1/2 

v(r)=-j-(I  -  I) 

'  m  I  '  a  c’ 


(123) 


which  when  substituted  into  the  relativistic  factor 
evaluated  at  the  sheath  edge  results  in: 


-r  =1/1 
m  a  c 


(124) 


So,  for  all  laminar  flows  t  =  I  /I  at  the  sheath  edge! 

mac 

3.4  Comparison  of  Constant  Density  and  Brillouin  Flow 
Theory  Predictions  to  MASK  Simulations 
In  many  short  cylindrical  MITL  simulations, 
non-uniformities  in  the  electric  fields  (usually  E^)  result 
in  laminar  flows  with  density  profiles  that  are  nearly 
constant.  In  the  following,  such  MASK  simulations  are 
compared  against  the  constant  density  and  Brillouin  flow 
theories.  The  magnetic  fields  for  three  cases  are 
presented:  1)  strongly  insulated  —  little  space  charge 

flow  (figure  12),  2)  fairly  well  insulated  —  more  space 
charge  flow  (figure  13),  and  3)  insulation  breaking  down  — 
greatest  amount  of  space  charge  flow  (figure  14). 

In  figure  12,  the  theoretical  results  essentially  lie 
on  top  of  the  MASK  results.  The  space  charge  current  flow 
(I  =  I  -  I  )  is  only  11%  of  the  cathode  current.  The 

6  H  C 

reason  the  Brillouin  and  constant  density  solutions  match  so 
closely  is  because  there  is  very  little  space  charge  current 


flow.  The  magnetic  field  is  primarily  due  to  the  cathode 
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current.  The  space  charge  current  distribution,  although 
radically  different  between  the  two  cases,  is  so  narrow  in  r 
that  the  total  current  (I  )  is  included  at  a  radius  of  .007 
meters.  By  Ampere's  law  the  magnetic  fields  must  match  from 
that  point  on.  Therefore,  the  fields  essentially  lie  on  top 
of  each  other  from  .007  meters  to  the  anode  (at  .05  meters) 
as  expec  ted . 

However,  in  figure  13.  it  is  clear  that  the  constant 
density  result  is  more  accurately  matching  the  simulation 
results  than  the  Brillouin  flow  theory.  It  is  also 
important  to  note  that  the  Brillouin  flow  solution  is  not 
satisfying  Ampere’s  law  at  the  cathode  or  anode.  The  reason 
for  this  is  that  the  Brillouin  flow  solution  is  valid  for 
only  certain  anode  and  cathode  currents  given  a  voltage 
across  the  diode.  These  currents  are  given  by  equations 
(61)  and  (62)  or  (101)-(103),  and  (124).  The  values  of  the 
anode  and  cathode  currents  in  figure  13  come  from  MASK  and 
do  not  satisfy  the  Brillouin  flow  restrictions.  A  nice 
feature  of  the  constant  density  theory  is  that  it  will 
produce  a  self  consistent  flow  solution  irregardless  of  the 
anode  and  cathode  currents  given. 

Figure  14  presents  the  magnetic  field  for  the  case  of  a 
semi-trapped  MITL  flow.  Once  again,  the  constant  density 
theory  corresponds  more  closely  to  the  MASK  calculation  than 
the  Brillouin  flow  solution.  The  match  between  the  constant 
density  theory  and  MASK  in  figure  14  is  not  as  good  as  it  is 
in  figure  13  because  the  density  profile  in  14  is  less 
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MASK  CALCULATION  M420  vs  THEORY 


Figure  13.  Theory  Vs.  MASK  Simulation  for  the  Case  of  a  Fairly 
Significant  Space  Charge  Flow. 
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constant  than  in  13.  One  should  also  note  that  the 

Brillouin  flow  solution  is  once  again  not  satisfying 

Ampere’s  law.  If  it  were,  it  would  match  the  constant 

density  solution  at  the  anode  and  the  cathode.  The  reason 

they  should  match  is  because  the  same  I  and  1  were  used  in 

a  c 

both  models.  The  reason  the  constant  density  solution  does 
not  exactly  match  MASK  at  the  anode  and  cathode  in  this  case 
is  because  the  MASK  magnetic  field  profile  seen  here 
corresponds  to  one  slice  across  the  flow  ( z=cons tan t ) .  For 
the  theoretical  treatments  the  cathode  and  anode  currents 
were  obtained  by  averaging  the  values  obtained  from  three 
slices. 

This  concludes  the  derivation  and  verification  of  the 
general  theory  of  laminar  flows.  In  the  next  chapter,  the 
general  theory  will  be  used  to  derive  an  electron  flow 
solution  downstream  of  a  PEGS.  However,  before  this  is 
done,  the  history  and  workings  of  a  PEGS  will  be  given. 
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CHAPTER  4 

PLASMA  EROSION  OPENING  SWITCH 

The  plasma  erosion  opening  switch  (PEGS)  has  been  used 
as  a  fast  output  switch  on  pulsed  power  generators.  It  can 
conduct  mega-Ampere  levels  of  current  for  lO's  of 
nanoseconds  before  switching  the  current  to  a  load  in  less 
than  10  nanoseconds.  The  PEOS  can  be  used  on  tera-Watt 
level  systems  for  inductive  storage,  pulse  compression,  rise 
time  sharpening,  and  power  multiplication  by  pulse 
compression  [27,31].  Indeed,  it  has  many  applications,  one 
of  which  was  mentioned  in  the  introduction  to  this 
dissertation  —  inertial  conf inement  fusion. 

In  cylindrical  geometries,  the  switch  typically 
conducts  current  through  a  plasma  injected  through  a  screen 
in  the  outer  anode  onto  a  central  cathode  (see  figure  15). 
When  the  ion  current  through  the  switch  exceeds  the  ion  flow 
supplied  by  the  injected  plasma  fill,  the  plasma  erodes  away 
from  the  cathode  and  the  switch  begins  to  open  [56].  As  it 
opens,  there  is  initially  a  significantly  higher  impedance 
across  the  downstream  power  feed  and  load  than  across  the 
switch.  This  mismatch  causes  electrons  to  continue 
streaming  from  the  cathode  to  the  eroding  PEOS  making  them 
unavailable  to  power  the  load.  As  the  switch  continues  to 
open  (the  injected  plasma  continues  eroding  away  from  the 
cathode)  the  impedance  mismatch  decreases  and  magnetic 
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insulation  insures  power  is  not  lost  across  the  switch 
itself. 

However,  between  the  switch  and  load,  electrons 
sometimes  continue  to  shunt  power  from  the  load  as  a  complex 
function  of  the  perturbing  PEGS  plasma,  applied  voltage, 
currents,  and  line  and  load  impedances.  This  potential  loss 
of  magnetically  insulated  electron  flow  is  a  major  concern, 
because  it  not  only  represents  a  power  loss,  but  it  can  also 
cause  the  PEGS  to  not  open  completely  (or  close 
prematurely).  Avoiding  the  regimes  in  which  losses  occur 
embodies  Chapter  5  and  modeling  the  electron  flow  in  the 
loss  free  region  embodies  part  of  this  chapter. 

Chapter  4  begins  with  a  historical  summary  of  PEGS 
research.  It  also  reviews  the  basic  physics  and  phases  of 
an  opening  PEGS.  Then,  the  electron  flow  downstream  of  a 
PEGS  during  the  opening  phase  is  studied.  In  regimes  of 
laminar  like  flows  downstream  of  the  PEGS,  the  geneial 
theory  (developed  in  Chapter  3)  is  used  to  derive  a  laminar 
flow  solution. 

4.1  Historical  Summary  and  Background  Physics 

This  is  a  historical  summary  of  PEGS  research  done  in 
the  United  States  [27].  Other  nations  such  as  the  Soviet 
Union  [26],  Japan  [35],  West  Germany  [36],  etc.  are  also 
involved  and  contributing  in  this  area,  but  to  a  lesser 
extent  at  present.  A  PEGS  was  first  used  in  1977  by  Mendel 
et.  al.  of  Sandia  for  prepulse  suppression  on  the  Proto  I 
Generator  [37].  In  1981,  Stringfield,  et.  al.  [3S]  used  a 
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PEGS  in  conjunction  with  experiments  done  on  the  Pithon 
Generator  at  Physics  International  Company.  Also  in  19S1 , 
the  Naval  Research  Labo ra to ry ( NRL )  became  involved  with  the 
PEGS.  Co  1 1 abo ra t i ve 1 y  with  C.  W.  Mendel  of  Sandiu,  they 
were  successful  in  suppressing  the  prepulse  and  preventing 
early  gap  closure  on  experiments  performed  on  the  Aurora 
Generator  at  Harry  Diamond  Laboratories  [39].  Ihe  NRL  group 
realized  at  this  point  that  a  PEGS  could  possibly  be  used  in 
an  inductive  storage  system  to  multiply  the  pulsed  power  by 
compressing  tfe  pulse  in  time. 

This  revolutionary  new  use  of  the  PEGS  was  demonstrated 
by  Meger,  et.  al.  of  NRL  in  1983  [dO].  Experiments  were 
performed  on  Gamble  I  of  NRL  that  showed  inductive  energy 
storage,  pulse  compression  and  power  multiplication.  Gamble 
I,  running  with  a  PEGS,  had  peak  powers  greater  tlian  the 
traditional  maximum  peak  powers  predicted  using  matched 
loads.  This  exciting  result  led  to  PEOSs  being  installed  on 
the  Super  Mite  Generator  at  Sandia  [dl]  and  on  PBFAII. 

Also,  high  current  experiments  were  performed  on  the  Black 
Jack  5  Generator  at  Maxwell  Laboratories  [d2].  All  these 
efforts  have  confirmed  whe  validity  of  t  li  i  s  new  use  for  the 
PEGS  . 

In  addition,  PEGS  work  is  now  being  done  by  NRL  in  the 
area  of  long  conduction  times  (.6  to  1.  it.  icrosecond)  [d3]. 

Ill  these  efforts,  the  switch  conducts  current  for  times  on 
the  order  of  a  mi  <':ro. second  before  swi  telling  the  current  to  a 


load  in  about  100  iianoseconds  . 


It  has  also  been  noted  that 
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in  pinch-reflex  ion  diode  experiments  performed  on  Gamble 

li.  ion  beam  brightness  is  substantially  increased  when  run 

with  a  PEGS.  This  is  attributed  to  faster  ion  turn-on 

There  are  many  other  PEGS  applications.  Gne  reported  on 

recently  at  the  1987  Plasma  Science  IEEE  Meeting  is  the  use 

of  PEGSs  for  repetitively  pulsed  power  applications  [45]. 

In  order  to  more  clearly  understand  the  possible 

applications  and  limitations  of  the  PEGS,  one  needs  to 

understand  the  basic  physics  and  phases  of  opening. 

The  phases  of  PEGS  opening  are  shown  in  figure  15. 

These  phases  and  the  physics  involved  are  presented  in  the 

remainder  of  this  section.  In  figure  15,  one  notes  four 

phases  of  PEGS  opening  which  in  accordance  to  NRL  are  [46]: 

a)  the  conduction  phase,  b)  the  erosion  phase,  c)  the 

enhanced  erosion  phase,  and  d)  the  magnetic  insulation 

phase.  A  brief  summary  of  each  phase  and  the  transition 

condition  from  phase  to  phase  is  given  next. 

During  the  conduction  pl.ase  the  resistance  across  the 

plasma  switch  is  approximately  zero.  A  thin  sheath  forms 

between  the  cathode  and  plasma.  Ion  current  from  the 

plasma  and  electron  current  (i  ,)  from  the  cathode  cross  the 

el 

sheath  in  a  bi-polar  space  charge  limited  fashion.  This  is 
given  by : 

1  /? 

i . / i  ,  =  { zm  /m . )  ,  (125) 

where  m.  is  the  ion  mass  and  m  is  the  electron  mass.  z  is 
1  e 

the  ion  charge  state.  The  electrons  that  are  emitted  from 
the  cathode  (in  a  SCL  fashion)  are  accelerated  from  the 


70 


cathode  towards  the  bottom  of  the  PEGS.  These  electrons 
enter  the  PEGS  plasma  with  a  large  amount  of  kinetic  energy. 
This  kinetic  energy  is  transformed  into  electromagnetic 
waves  that  cause  a  large  number  of  anomalous  collisions  or 
mixing.  The  end  result  of  such  mixing  could  be  anamolous 
magnetic  diffusion  and  wide  current  channels.  Also,  the 
current  would  be  predominately  radial  in  nature.  When 
anamoulous  collisions  are  not  taken  into  account  fluid  code 
simulations  in  1  and  2  dimensions  [^T],  PIC  code  simulations 
[48.  49]  and  theoretical  treatments  [50]  show  non-radial 

narrow  current  channels.  However,  experiments  done  at  NRL 
clearly  show  wide  current  channels  and  predominantly  radial 
current  paths  [51],  When  anamoulous  collisions  are  taken 
into  account  —  with  a  collision  frequency  on  the  order  of 
the  plasma  electron  frequency  —  PIC  [52]  and  fluid  [53,54] 
code  simulations  match  the  experiments  done  at  NRI,.  An 
instability  responsible  for  such  anamoulous  collisions  has 
been  suggested  and  analyzed  by  Kulsrud,  Ottinger  and 
Grossman  [55]  . 

In  general  m.  >>  m  .  so  that  the  current  crossing  the 
1  c 

switch  (i^  =  is  approximately  equal  to  i^j  (by 

(125)).  When  i^  becomes  so  large  that  i^  is  greater  than 

that  supplied  by  the  injected  plasma,  the  conduction  phase 

ends  and  the  erosion  phase  begins.  This  happens  when  [56]-' 
rm  z  1  1  /2 


^  e 

i  m  A 
P 

where  v,  is  the 
d 

numbe  r  density, 


i  >2TrIrnezv.  ,  (126) 

s  c  p  a 

plasma  drift  velocity,  n^  is  the  plasma 
1  is  the  axial  length  of  the  plasma  and  A  is 
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the  atomic  number. 

When  more  ions  are  pulled  off  the  bottom  plasma  surface 
than  the  plasma  drift  can  resupply,  the  plasma  erodes  away 
from  the  cathode.  This  is  shown  in  figure  15b.  The  growth 
of  the  gap  (D)  between  the  plasma  and  cathode  (figure  15b) 


i s  given  by  [56] : 


dP 
d  t 


i  .  -  n  zev  ,2Trr  1 
p  d  c 

n  ze27rr  1 


(127) 


P  c 

As  the  gap  grows,  current  begins  flowing  through  the 
load.  This  leads  to  partial  magnetic  insulation  of  the 
space  charge  electron  flow.  The  magnetically  trapped 
electrons  electrostatically  attract  ions  enhancitig  the 
plasma  erosion.  This  is  called  the  enhanced  erosion  phase 
in  which  the  increased  ion  current  is  given  by  [56, .57]: 


1  /9 

i.  =  i^(2zm^/m.)  (-r+ni/D 


(128) 


As  more  current  passes  through  the  load,  the  magnetic 
insulation  is  increased.  When  the  current  through  the  load 
is  so  large  that  no  r  rons  reach  the  anode,  the  magnetic 

insulation  phase  begins.  The  minimum  required  insulating 
cu-rent  for  this  as  proposed  bj'^  Ottinger,  et.  al.  is  [56]: 


i 

c  r  1  t 


1 . 6  X  8500 


l2L 


1  n  {  r 


1/2 

T 


(129) 


This  is  the  standard  relation  for  the  critical  current 
multiplied  by  a  factor  of  1.6.  The  factor  is  obtained  in 
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the  analysis  of  a  couple  calculations  done  in  the  2MV  regime 
[58]  in  conjunction  with  the  experimentally  measured 
impedance  behaviour  of  a  pinch  reflex  diode  by  NRL .  A  more 
in  depth  numerical  and  theoretical  analysis  is  performed  in 
Chapter  5  to  more  accurately  determine 
■4.2  Magnetic  Insulation  Phase  Electron  Flow 

When  the  current  through  the  load  exceeds  the  critical 
current  (Chapter  5),  the  electron  flow  is  insulated  away 
from  the  anode  and  the  magnetic  insulation  phase  is 
underway.  Even  though  the  flow  is  insulated,  it  can  be 
quite  chaotic  downstream  of  the  PEGS.  In  such  a  situation 
there  are  various  instabilities  [11,22,32.59-63]  that  can 
complicate  the  understanding  of  the  electron  dynamics. 
However,  in  the  MASK  calculations  done  in  this  research  they 
played  a  negligible  role.  This  could  be  because  the  length 
of  the  MITL  between  the  switch  and  load  is  so  short  that  the 
instabilities  do  not  have  a  chance  to  develop. 

An  interesting  aspect  of  having  a  short  MITL  between 
the  switch  and  load,  is  that  the  magnetic  insulation  of  the 
system  can  be  made  to  occur  on  a  faster  scale.  This  can  be 
understood  by  considering  an  electromagnetic  pulse 
propagating  down  an  unperturbed  MITL.  At  the  pulse  front, 
electrons  flow  directly  from  the  cathode  to  the  anode. 
Downstream  of  the  pulse  head,  electrons  are  bent  into  orbits 
that  are  governed  by  the  self-limited  impedance  of  the  line 
(figure  2a).  This  state  continues  until  the  pulse  head 
strikes  the  load.  If  the  impedance  of  the  load  verses  the 
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impedance  of  the  line  as  a  function  of  voltage  is  less  than 
the  ratio  given  in  figure  39,  then  an  insulation  wave  will 
propagate  from  the  load  end  of  the  system  back  down  the 
line.  This  insulation  wave  causes  an  increased  magnetic 
field  that  more  tightly  traps  the  electron  space  charge  flow 
to  the  cathode  (see  figure  16).  The  longer  the  line  or 
distance  between  the  load  and  PEGS  the  longer  it  will  take 
for  this  increased  insulation  of  electron  flow  to  take 
effect.  If  the  load  impedance  changes  as  a  function  of  time 
the  effect  on  the  PEGS  will  be  retarded  proportional  to  its 
distance  from  the  load. 

It  is  interesting  to  note  that  in  PEGS  experiments  done 
at  NRL,  it  has  been  found  that  minimizing  the  distance 
between  the  PEGS  and  load  reduces  current  losses  and  results 
in  more  efficient  switching  [64].  This  is  probably 
attributable  to  reduced  instabilities  and  quicker  magnetic 
insulation  in  shorter  systems. 

Before  one  considers  the  instabilities  of  the  insulated 
flow  downstream  of  the  PEGS,  the  fundamentals  of  the  flow 
should  be  understood.  In  MASK  simulations  of  a  stationary 
PEGS,  the  fundamental  properties  of  the  flow  have  been 
examined.  The  PEGS  perturbation  launches  the  electron  flow 
in  a  manner  consistent  with  an  drift.  The  electron 

orbits  are  often  laminar  in  nature.  However,  for  well 
trapped  magnetically  insulated  flows  downstream  of  a  PEGS, 
the  laminar  flow  introduced  by  the  perturbations  can  bo 
overshadowed  by  the  magnetic  fields.  When  this  happens,  the 
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Figure  16.  Insulation  Wave  in  a  Load-Limited  Line. 

The  insulation  wave  propagates  from  the  load  end 
back  up  the  line  resulting  in  e  more  tightly 
trapped  electron  flow. 


flow  enters  non-laminar  electron  orbits  that  are  typically 
prolate  or  curtat,.  In  this  section,  the  launching 
mechanism  and  examples  of  laminar  and  non-laminar  flows 
downstream  of  an  abruptly  varying  structure  {with  and 
without  ion  emission)  are  presented. 

4.2.1  Beam  Launching 

In  MASK  simulations,  the  PEGS  often  launches  the 
electron  flow  (from  the  cathode  and)  into  the  cathode/anode 
gap.  An  example  of  this  is  shown  in  figure  17  which  is  a 
vector  plot  of  the  current  densities  for  problem  PN06  (Table 
C4  of  Appendix  C) .  In  simulation  PN06,  the  load  is  fairly 
large  and  results  in  a  poorly  trapped  flow  downstream  of  the 
PEGS. 

In  detailed  analysis  of  various  MASK  simulations,  it 
appears  that  the  primary  mechanism  for  launching  the 
electrons  is  a  component  of  the  ^x?  drift  in  the  positive 
radial  direction.  The  ?x^  drift  velocity  is  given  by: 

/•  1  ^  ^  2  ,  ?x]?  f  1  rirw 

V  =  (1  +  — V  )  — ^  (130) 

e  a  gz: 

Since  ^  =  -B(r)6.  ?  must  have  a  significant  component  in  z 
for  beam  launching.  Upstream  of  the  structure  (PEGS)  there 
is  a  -z  component  of  ?  which  suppresses  launcliing;  however, 
downstream  of  the  PEGS  ^  has  a  +z  component.  One  should 
also  note  that  in  going  from  a  region  of  no  space  charge 
flow  (i.e.  an  insulator)  to  one  with  space  charge  flow  (see 
figure  1)  there  will  also  be  a  +z  component  of  ?  that  will 
cause  partial  or  complete  beam  launching.  Complete  beam 
launching  refers  to  launched  electron  flows  that  completely 


leave  the  cathode. 


To  obtain  a  more  quantitative  idea  of  this,  one  may 
electrostatically  determine  the  fields  for  the  structure  of 
interest.  The  geometry  to  be  considered  is  given  in  figure 
18.  To  simplify  the  mathematics  and  conserve  much  of  the 
physics  one  can  consider  the  region  within  the  dashed  lines 
in  cartesian  coordinates.  Physically  this  is  equivalent  to 
letting  r  ,  r  go  to  infinity  while  maintaining  d  and  D. 

3i  C 

The  expressions  obtained  from  an  analysis  similar  to  the 
following,  but  done  in  cylindrical  coordinates,  are 
cumbersome . 

In  cartesian  coordinates,  the  potential  { )  at  y  =  D  is 
V  ,  a t  y  =  0  i t  i s  0 ,  and  at  infinity  the  potential  is 
bounded.  Atx  =  0,  0<y<d,  if  =  V  y/d .  By  standard 

Si 

techniques  it  is  found  that  the  potential  within  the  dashed 
region  is: 
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where  ^  .  The  negative  of  the  gradient  of  (131)  gives 

the  electric  field  components: 

sin(A^y )cxp(-X^x) 

(132) 
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Figure  18.  Geometry  for  Theoretical  Analysis  of  Beam  Launching. 
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The  magnetic  induction  can  be  approximated  using  Ampere's 


law  -  given  the  current  through  the  load  (1  ): 

pi 

g  ^  . 

w 

where  w  is  the  width  of  the  system. 

From  (130)  becomes: 
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,nd 


Since  the  finite  Larmour  radius  effect  is  often  a  2  order 
effect,  it  is  instructive  to  write  out  the  drift  velocity 
components  neglecting  the  terms  related  to  the  finite  radius 
effect.  From  (132)-(135): 
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,  n+  1 


sin(X^y)cxp(-X^x) 
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In  the  limit  x 
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"  .  (136)  and  (137)  go  to: 
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z  '  z  ' 

Vy  =  0  .  (139) 

This  is  as  expected.  Far  away  from  the  PEOS  perturbation, 
the  electrons  ^x^  drift  only  in  the  x  direction  (figure  18). 

However,  close  to  the  perturbation,  there  is  a  sizable 
+  y  component  of  drift  (see  equation  (137)).  Tliis  y 
component  of  drift  has  a  maximum  close  to  the  perturbation 
edge.  It  decreases  to  zero  at  the  anode  and  the  catliode. 

It  also  decreases  exponentially  as  one  moves  away  from  the 

As 

perturbation  in  x.  Therefore,  after  the  perturbation 
launches  the  beam,  it  has  less  and  less  effect  on  the 
receeding  electron  flow.  In  addition,  since  the  y  component 
of  drift  goes  to  zero  at  the  cathode  an  initially  wider 
space  charge  electron  flow  (  not  as  well  trappd  )  will  be 
more  severely  effected  (Note:  x-»z  .  y-»r  ,  z->0  when 

approximating  cylindrical  geometries). 

To  more  realistically  see  how  this  works,  a  MASK 
calculation  is  presented  in  which  the  ?x^  launching 
mechanism  is  studied  for  an  insulated  system  with  an  ion 
emitting  structure  (similar  to  an  opening  PEOS).  The  study 
requires  one  to  first  estimate  the  average  components  of  the 
and  ^  fields  across  the  flow  in  the  launching  region. 

From  figure  19,  the  average  magnetic  induction  across  and 
along  the  electron  stream  in  the  launch  region  is  (the 
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launched  electrons  are  within  the  dashed  lines  of  figure 
19): 

^  S  -1.8  e  (Tesla) . 

From  figure  20,  the  average  r  component  of  the  electric 
field  along  the  flow  in  the  launch  region  is  approximately: 
S  -3.5x10®  r  (V/m) . 

From  figure  21,  the  average  z  component  of  the  electric 
field  is  seen  to  be  approximately: 

=  7.8x10^  2  (V/m) 

Equation  (130)  implies  that  the  ^x^  components  of  drift 
ve 1 oc i ty  ar e : 


V 

r 


(  m/sec ) . 


V  =  ^  S  1.9x10®  (m/sec) . 

Z  Kg 

From  figures  22a  and  22b  it  is  clear  that  the  launching 

of  the  beam  occurs  primarily  between  z  =  .07  and  z  =  .11  (Az 

=  .04  Meters).  Using  the  calculated  velocities,  it  is 

possible  to  estimate  the  sheath  d i sp 1 acemcn  i  ( A s ) . 

A  7 

As  =  —  V  =  .009  Meters 
V  r 


z 

This  displacement  is  approximately  equal  to  that  predicted 
by  MASK  as  shown  in  figure  22a . b .  Since  one  can  predict  the 
sheath  displacement  by  assuming  an  ?xB  drift,  this  implies 
the  ?x^  drift  is  the  primary  launching  mechanism.  Other 
calculations  were  examined  with  and  without  ion  emission 
from  their  respective  structures.  In  every  case  the 
electron  trajectories  near  the  perturbation  are  consistent 
with  an  ?x?  launching  mechanism. 


f/t  20/3 
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Figure  19.  Magnetic  Field  V.’ithin  the  Launched  Beam. 

The  average  magnetic  field  on  the  electrons  at  the  right 
edge  of  the  PEGS  is  approximately  given  by  contour  level 
number  6.  The  electrons  move  between  the  dashed  lines 
in  the  general  direction  of  the  arrow. 


P3 


2  -5.9x10®(V/i7i) 

3  -4.8x10® 

4  -3.7x10® 

5  -2.56x10® 

6  -1.44x10® 

7  -3.33x10^ 


Figure  20.  Electric  Field  (r)  Within  the  Launched  Beam. 

The  approximate  r-component  of  the  electric  field 
acting  on  the  launched  electrons  (between  the  dotted 
lines)  is  given  by  averaging  contour  levels  3,  4 
and  5 . 
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Figure  21.  Flectric  Field(z)  Within  the  Launched  Beam. 

The  average  z-component  of  the  electric  field  in  the 
launch  region  is  given  approximately  by  contoi.r  level 
number  8  (it  bisects  the  electron  flow) .  The  electrons 
flow  between  the  dashed  lines  in  the  direction  indicated 
bv  the  arrow. 
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■1.2.2  Laminar  Flows 

If  the  impedance  of  the  load  relative  to  the  impedance 
of  the  line  as  a  function  of  voltage  is  within  approximately 
10%  of  the  ratio  give  in  figure  39.  the  launched  electron 
beam  will  propagate  in  a  laminar  fashion  (see  Appendix  C, 
Tables  C1-C9).  This  of  course  applies  only  to  the  electron 
flow  downstream  of  an  abruptly  varying  structure  as  sliown  in 
figure  23.  The  flow  referred  to  as  being  laminar  here  is 
not  really  laminar  in  that  f  ^  0.  However,  the  electrons 
typically  travel  10  or  more  Larmour  radii  in  the  z  direction 
per  gyration.  Since  the  excursion  in  r  is  close  to  two 
Larmour  radii  per  gyration  the  orbits  are  fairly  laminar. 

In  the  following  these  semi-laminar  orbits  will  be 
investigated  for  one  MASK  run  (characteristic  of  many). 

Then  using  the  general  laminar  flow  theory,  a  flow  solution 
is  generated  for  laminar  electron  flows  downstream  of  the 
PEGS.  The  laminar  flows  are  of  particular  interest  since  it 
is  often  desirable  to  run  systems  with  as  large  a  load  as 
possible  and  still  maintain  magnetic  insulation  (this  puts 
one  into  the  laminar  regime). 

The  simulation  to  be  considered  is  MASK  run  PN06  (see 
Appendix  C,  Table  C5).  The  maximum  potential  difference 
between  the  anode  and  cathode  is  2MV.  There  is  no  ion 
emission  from  the  perturbation  and  the  problem  has  been  run 
to  equilibrium.  In  figure  23.  one  can  see  the  typical 
semi  -  laminar  electron  orbits  under  consideration  (Note  that 
if  the  r  -ind  z  axis  nad  the  same  scale  the  plot  would  appear 
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4  times  longer  in  the  z  direction).  From  figure  23.  it  can 
be  seen  that  the  electrons  move  5  to  6  times  further  in  the 
z  direction  than  their  excursions  in  the  r  direction  per 
gyration.  The  orbits  are  fairly  laminar.  This  is  due  to 
the  perturbation  modifying  the  electron  orbits  and 
velocities  such  that  the  electric  and  magnetic  field  forces 
cance  1 

=  0.  .  (140) 

If  E  <  V  B„  the  electrons  are  accelerated  towards  the 
r  z  0 

cathode;  however,  if  E  >  v  the  electrons  arc  accelerated 

r  z  0 

towards  the  anode.  If  E  =  v  the  electrons  move  in 

r  z  6 

perfectly  laminar  orbits.  A  relative  measure  of  the  laminar 

nature  of  the  electron  flows  can  be  obtained  by  comparing 

the  values  of  E^-v^Bg  at  the  turning  point  for  different 

load  impedances.  The  more  laminar  flows  will  have  smaller 

absolute  values  for  E  -v  B„.  This  will  be  demonstrated  by 

r  z  G 

examining  two  MASK  simulations. 

The  outer  turning  point  for  the  electrons  in  MASK 
simulation  PN06  is  approximately  given  by  a  radius  of  .045 

g 

meters  {r  ).  From  figures  24  and  25,  E  (r  )  =  -1.6x10  and 
'  m  ^  ^  r  m ' 

5  '’l 

)  =  “1^  (5.0x10  ).  From  figure  26.  v  (r  )  =  —  = 

0'm''  ^o'  '  z  m^  1 

2.8x10^.  This  implies: 

E  -  V  Bq  =  1 .67x10^  .  (  141  ) 

r  z  o 

Multiplying  this  by  the  electron  charge  results  in  a 
negative  r  directed  force  which  causes  the  electrons  to  move 
towards  the  cathode.  However,  it  is  a  relatively  small 
force  as  can  be  seen  in  comparison  to  the  next  simulation. 
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The  next  case  is  for  simulation  PN04  (Appcr  dix  )  . 

This  case  corresponds  to  simulation  PN06  in  everyway  except 

that  the  load  impedance  is  reduced  which  causc's  a  more 

tightly  trapped  non-laminar  electron  flow  downsiream  of  the 

perturbation.  The  typical  electron  orbits  for  simulation 

PN04  are  shown  in  figure  27.  The  outermost  point  for  the 

electron  orbits  is  r  =  .045  Meters.  From  figuies  2S  and  29 

m 


the  fields  are: 
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E  { r  )  S  -1 . 375x10 
r  m ' 

Be(r,.)  =  -.82 

From  figure  30,  the  z  component  of  velocity  is  v 


E  -  V  B.  =  +1  .xlO 
r  z  6 


2.87x10  .  So. 

(142) 

In  comparing  equation  (141)  and  (142),  the  minus  r 
directed  force  of  (142)  is  six  times  greater  than  that  given 
in  equation  (141).  Considering  this  result,  it  is  not 
surprising  that  the  electron  flow  downstream  of  the  PEOS  in 
simulation  PN06  is  laminar  relative  to  PN04. 

As  mentioned  previously,  it  is  often  desirable  to  power 
loads  with  the  largest  possible  impedance,  resulting  in 
semi-laminar  flows  downstream  of  abrupt  perturbations  (i.e. 
PEOS).  From  a  series  of  PEOS  calculations  (Tables  C7-C9)  it 
is  found  that  the  electron  density  profiles  downstream  of 
the  perturbation  are  often  approximately  given  as  follows: 


p( r )  =  0 

p ( r )  =  cons  tan  t 
p(r)  =  0 


r  <  r  <  r  . 
c  “  ~  1 


r  .  <  r  <  r 
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r  <  r  <  r 
o  "  ~  a 


(143) 
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where  r  is  the  cathode  radius,  r.  is  the  inner  radius  of 
c  1 

the  launched  be:.ii.  r  is  the  outer  radius  of  the  beam  and. 

o 

r^  is  the  anode  radius.  Since  the  flow  is  primaiily  laminar 
the  general  theory  of  laminar  flows  (from  Chapter  3)  can  be 
used  to  predict  flow  quantities  of  interest.  Since  the 
density  is  approximately  constant  across  the  flow,  the 
constant  density  solution  can  be  modified  to  give  the 
downstream  electron  flow. 

The  laminar  flow  solution  downstream  of  a  PEOS  is  given 

in  equations  ( 1 43 ) - ( 1 62 ) : 

Forr  <r<r.  , 

c  ~  -  1 

=  t(r)  =  <p(r)  =  t{r)  =  f(r)  =  0  ,  (144) 
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For  ^  ^  ^  (see  equations  ( 1 05 ) - ( 1 1 4  )  )  : 
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where  a  =  c_/4  ,  b  =  -c_r.  .  c,  =  c_r.  +  1 
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For  r  <  r  <  r  (see  equations  ( 1 1 5 ) - ( 1 20 ) ) : 

O  3. 
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and 


v(r)  =  p(r)  =  7{r)  =  0  . 


where  ^  2  ' 

2r  . 
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The  electron  density  is  given  by 
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The  radius  to  the  outer  edge  of  the  beam  is  obtained 
itteratively  from  (equations  (155)  at  the  anode  and  (161)): 
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One  may  obtain  a  detailed  flow  solution  for  the  laminar  flow 
downstream  of  the  PEGS  from  these  equations  given  the 
applied  voltage  across  the  diode  (V  ).  the  anode  (I  )  and 

^  Cl 

cathode  (I  )  currents,  the  anode  (r  )  and  cathode  (r  ) 

C  3.  C 

radii  and  the  inner  radius  (r.)  of  the  launched  beam. 

All  these  variables  are  easily  obtained  except  the 
inner  radius  of  the  beam.  The  inner  radius  (t^)  depends  on 
the  size  of  the  perturbation,  diode  voltage,  anode  and 
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cathode  currents,  and  geometry.  How  these  relate  to  the 
inner  beam  radius  has  not  been  quan titatively  dctertnined. 

Oiie  can  estimate  r.  based  on  the  expected  electric  fields  in 
the  launch  region  as  was  done  in  section  4.2.1.  The  inner 
beam  radius  for  the  launched  beam  can  be  very  roughly 


approximated  by: 


4.2.3  Non-Latninar  Flows 


(163) 


It  is  important  to  stress  that  if  the  impedance  of  the 
load  relative  to  the  impedance  of  the  line  as  a  function  of 
voltage  is  less  than  that  which  gives  laminar  flows  (see 
section  4.2.2),  the  launched  beam  will  form  circulating 
rings  of  charge  (figure  31a. b).  These  rings  will  ^x^  drift 
down  the  transmission  line.  The  individual  electron  orbits 


are  usually  prolate  and  curtate  (see  figures  1  and  31b). 

In  highly  trapped  electron  flows,  the  launched  beam 
does  not  come  close  to  satisfying  (140)  Since  the  electron 
beam  curls  up,  it  no  longer  can  completely  shield  the 
cathode  from  electric  fields.  Therefore,  between  vortex 
rings  (circulating  electrons,  ?x?  drifting  together  down  the 
line)  space  charge  limited  electron  emission  occurs  from  the 
cathode.  Underneath  each  vortex  ring,  there  is  no  electron 
emission  because  of  the  rings  shielding  effect  (see  fig 
33a).  This  phenomena  is  rather  irregular  and  can  be  most 
easily  predicted  using  a  2-dimensional  electromagnetic  code. 
The  general  theory  of  laminar  flows  can  not  be  used  in  such 
a  si tua t i on . 
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Figure  27.  Non-Laminar  Electron  Orbits 
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CHAPTER  5 


MAGNETIC  INSULATION  THRESHOLD 
There  are  three  types  of  magnetically  insulated  flow; 

1  )  1 oad- 1  imited.  2)  self-limited,  and  3)  cons  tan t- f 1 ux 
-limited  flow  (see  figure  2).  In  load  limited  cylindrical 
MITLs  (see  i  i^xyte  2).  there  is  a  maximum  load  impedance  (as 
a  function  of  voltage)  above  which  magnetic  insulation 
degenerates  across  the  load.  This  results  in  current  being 
excluded  from  the  load  as  shown  in  figure  32.  This  loss  of 
current  is  undesirable.  To  determine  what  the  magnetic 
insulation  threshold  is  (above  which  current  is  lost)  for  a 


of  this  Chapter.  —  r— r 

Now.  it  is  interesting  to  note  that  a  load  limited  MITL 
transitions  to  a  self-limited  MITL  at  the  point  current 
begins  to  be  excluded  from  the  load.  This  is  fortunate 
because  current  flows  in  self-limited  magnetically  insulated 

4 

y 

lines  httve  been  treated  theoretically  in  cylindrical  ^ 
coordinates  by  Wang  and  DiCapua  [33]  assuming  a  Brillouin 
flow  field.  Wang  and  DiCapua  determine  the  minimum  energy 
operating  point  of  self-limited  MITLs  as  a  function  of 
applied  voltage.  Realizing  that  the  normalized  enefrgy  of 
the  Brillouin  flow  at  the  transition  point  (between 
load-limited  and  seif-limited  flows)  is  inversely 
proportional  to  the  normalized  load  impedance  allows  one  to 
invert  Wang  and  DiCapuas’  results  to  obtain  figure  33.  So. 


ADDENDUM:  FJSRL-TR-89-0008,  Aug  89,  Magnetically  Insulated  Electron  Flows  in  Pulsed  Power 
Systems,  by  R.I.  Lawconnell  103 
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Figure  32.  Current  Loss  Due  to  Magnetic  Insulation  Deterioration. 
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Figure 33.  Magnetic  Insulation  Threshold  for  Brillouin  Flows. 

This  is  an  inversion  of  Wang  and  DiCapuas'  minimum 
energy  operating  point  verses  voltage  graph  (33). 
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instead  of  an  energy  verses  voltage  plot,  a  normalized  load 

i  tnpedance  ( z  j /z^ )  verses  applied  voltage  plot  is  obtained, 

where  z,  and  z  are  the  load  and  line  impedances 
1  o 

respectively.  In  this  work,  such  a  plot  is  referred  to  as 
the  magnetic  insulation  threshold  of  the  system. . 

In  the  following,  it  becomes  apparent  that  all  load 
limited  lines  with  trapped  flows  operate  below  the  Zj/z^ 
magnetic  insulation  threshold  curve.  MITLs  operating  above 
the  curve  shown  in  figure  33  experience  current  loss  (as 
shown  in  figure  32).  Thus,  figure  33  becomes  a  guide  to  the 
magnetic  insulation  threshold  for  infinitely  long  (no 
significant  end  effects  observed  in  the  flow)  cylindrical 
MITLs  with  no  perturbations.  This  curve  is  justified  by 
detailed  MASK  calculations  in  this  chapter.  In  addition, 
curves  similar  to  that  of  figure  33  are  computationally 
generated  in  which  the  magnetically  insulated  tlireshold  for 
finite  systems  with  structures  (some  structures  simulate  an 
opening  PEGS)  are  presented.  Some  of  these  magnetic 
insulation  threshold  curves  are  justified  by  using  them  to 
successfully  predict  experimental  results  from  a  Gamble  II 
PEGS  experiment. 

5.1  Slightly  Perturbed  Flow 

In  this  portion  of  the  magnetic  insulation  threshold 
study,  data  was  pulled  from  approximately  fourty  MASK 
simulations  (Appendix  C,  Tables  C1-C3)  of  finite  MITLs  (see 
figure  3a).  This  data  is  plotted  and  discussed  in  the 
following. 

The  simulations  typically  require  about  4  nanoseconds 
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of  physical  time  for  the  MITL  to  come  to  equilibrium.  All 
the  values  in  Tables  C1-C3  correspond  to  the  equilibrium 
condition.  Longer  MITLs  would  require  more  time  to  come  to 
equilibrium  as  discussed  in  section  4.2  and  indicated  on 
f i gure  1 6 . 

The  threshold  for  complete  magnetic  irisulaiion  is  best 
represented  by  a  Zj/z^  vs  voltage  graph,  where  z^  is  the 
load  impedance  and  z^  is  the  transmission  line  vacuum 
impedance.  The  voltage  is  that  applied  across  the 
transmission  line.  A  typical  data  point  for  such  a  graph  is 
obtained  by  doing  3  or  4  short  MASK  calculations  in  which 
the  diode  gap  (figure  3a)  is  increased  while  maintaining  a 
constant  voltage  difference  across  the  line.  When  the 
current  begins  to  be  excluded  from  the  load  (figure  32), 
magnetic  insulation  of  the  load  deteriorates. 

Determining  when  this  begins  is  easily  done  if  one 

plots  Zj/z^  vs  the  diode  gap  spacing  (d^).  This  results  in 

a  line  of  constant  positive  slope  as  long  as  the  flow  is 

trapped.  However,  when  d^  becomes  so  large  that  current 

begins  to  be  excluded  from  the  load,  the  slope  of  z,/z  vs 

1  o 

d  goes  to  zero  (where  z,  is  defined  to  be  V  /I  ).  An 
example  of  this  type  of  graph  is  given  in  figure  34.  The 
data  points  come  from  simulations  M407-H424  of  Table  Cl. 

The  point  at  which  the  slope  of  the  line  goes  to  zero  (in 
figure  34)  corresponds  to  the  magnetic  insulation  threshold 
for  the  given  conditions.  The  slope  goes  to  zero  because 
when  the  load  impedance  exceeds  the  self-limited  impedance 
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Figure  34.  Normalized  Load  Impedance  Vs.  Diode  Gap  Spacing. 

The  point  at  which  the  slope  of  the  line  goes  to  zero 

corresponds  to  the  magnetic  insulation  threshold.  This 

is  for  the  case  of  V  =2x10^  volts,  r_=.05,  rf.=  .025  meters. 

o  a 
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of  the  line,  the  current  flowing  through  the  line  begins  to 
short  across  the  line.  This  means  that  as  one  continues  to 
increase  the  diode  gap  spacing  more  and  more  current  is 
excluded  from  the  load.  However,  the  current  flowing  in  the 
system  (including  the  short)  is  given  by  the  self-limited 
impedance  of  the  line  —  which  is  a  constant. 

After  plotting  the  data  in  Tables  C1-C3,  figure  35 
results.  It  should  be  noted  that  each  point  on  this  graph 
is  found  by  the  technique  illustrated  in  figure  3^1  and 
discussed  above.  Figure  35  gives  the  magnetic  insulation 
threshold  in  terms  of  z^/z^  vs  voltage  plots.  The  upper 
line  corresponds  to  coaxial  MITLs  with  no  perturbations.  It 
indicates  that  transmission  lines  opperated  below  the  line 
lose  no  current  out  of  the  load.  The  triangled  dashed  line 
corresponds  to  MITLs  run  with  space  charge  limited  ion 
emission  from  the  anode.  Such  lines  must  be  run  at  a  lower 
load  impedance,  as  compared  to  the  ideal  case,  for  magnetic 
insulation  to  be  effective. 

It  is  also  clear  in  figure  35  that  varying  the 
anode/cathode  radii  had  an  insignificant  effect  on  the 
magnetic  inslulation  curves.  This  is  consistent  with  the 
flow  theory  of  Wang  [33]  mentioned  earlier.  In  addition, 
the  Brillouin  flow  self-limited  impedance  curves  are 
slightly  higher  than  the  ca 1 cu 1  a t i ona 1  curves;  but,  very 
close  to  them.  Since  all  orbits  are  possible  in  the  MASK 
calculations  and  the  flows  are  predominantly  non-B r i 1 1 ou i n 
(as  demonstrated  in  Chapter  2),  it  is  surprising  to  see  how 
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Figure  35.  Magnetic  Insulation  Thresholds  for  Systems  Without 
Perturbing  Structures. 

Cylindrical  transmission  lines  with  insignificant  ion 
emission  off  the  anode  lose  no  current  through  the 
load  if  operated  below  the  solid  line.  MITLs  with  space 
charge  limited  ion  flow  off  the  anode  lose  no  current 
if  operated  under  the  dotted  line. 


well  the  Brillouin  minimum  energy  theory  comes  to  predicting 
the  magnetic  insulation  threshold  found  in  the  ca 1 cu 1  a t i ona 1 
study  for  non-pe r tu r bed  MITLs.  The  magnetic  insulation 
thresholds  for  systems  with  structures,  with  and  without  ion 
emission  are  presented  in  the  following  section. 

5.2  Disrupted  Flows 

The  data  for  the  following  study  comes  from 
approximately  sixty  MASK  simulations  (Appendix  C.  Tables 
C4-C9).  The  MASK  mesh  for  these  simulations  is  given  in 
figure  3b.  All  the  calculations  were  run  to  equilibrium. 
Part  of  these  runs  (Tables  04-06)  were  done  with  a  large 
conducting  structure  in  the  mesh  that  did  not  emit  ions. 

The  remainder  of  the  runs  (Tables  07-09)  had  ion  emission 
from  the  structure  and  represented  PEOSs  in  various  stages 
of  opening.  As  a  result,  the  magnetic  insulation  thresholds 
that  are  obtained  apply  to  extremely  perturbed  systems. 

In  these  simulations  the  axial  extent  and  location  of 
the  perturbation  has  a  small  effect  on  the  magnetic 
insulation  properties  of  the  system.  The  vacuum  impedance 
of  the  line  in  the  switch  (perturbation)  region  is  the  most 
important  factor  governing  insulation.  In  fact,  the  vacuum 
impedance  through  the  switch  governs  the  insulation  of  the 
system  in  a  very  similar  fashion  as  the  vacuum  impedance  of 
the  MITL  in  the  last  section.  Because  of  this,  the  magnetic 
insulation  threshold  is  most  appropriately  represented  with 


a  z,/z  vs  voltage  plot,  where  z  is  the  vacuum  line 

Ip  &  e  .  p 

impedance  under  the  structure  (601n( r^/r^) ) . 

Once  again,  the  transition  from  trapped  to  untrapped 


flow  can  be  seen  in  the  Zj/z^  vs  plots.  A  representative 

Zj/Zp  vs  plot  is  given  in  figure  36  for  a  non-ion 

emitting  structure.  Figure  37.  is  a  compilation  of  z,/z  vs 

1  P 

voltage  data.  The  top  line,  connecting  diamonds  and 
squares,  corresponds  to  the  magnetic  insulation  tlireshold 
for  non-ion  emitting  structures. 

The  more  interesting  and  crucial  set  of  runs  for 
understanding  a  PEGS  are  represented  by  the  triangles, 
circles  and  crosses  connected  by  the  dashed  line.  In  these 
runs,  the  PEGS  structure  is  a  space  charge  limited  ion 
emitter.  If  the  pulsed  power  system  is  run  witli  a  load  and 
at  a  voltage  corresponding  to  the  region  below  the  dashed 
line  all  current  passes  through  the  load.  If  the  load  and 
voltages  are  such  that  the  system  operates  above  the  dashed 
line,  the  system  will  not  be  fully  magnetically  insulated 
and  current  will  be  lost.  In  the  next  section,  functions 
are  generated  that  fit  the  data.  These  functions  are  then 
used  in  section  5.4  to  model  a  Gamble  II  experiment. 

5.3  Fits  to  Magnetic  Insulation  Data 

It  is  desirable  to  have  a  continuous  function  that  fits 
the  data  linked  by  the  dashed  line  of  figure  37.  To  do 
this,  it  is  necessary  in  this  study  to  fit  the  data  for  the 
magnetic  insulation  of  a  transmission  line  with  no 
perturbation  (see  figure  35).  This  is  because  more  data 
points  are  available  for  the  unperturbed  case.  Since  the 
graphs  have  essentially  the  same  form,  the  unperturbed  case 
can  be  used  to  give  an  estimate  of  the  perturbed  case  in  the 
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Figure  37.  Magnetic  Insulation  Threshold  Data  for  Systems  with 
Perturbing  Structures. 

The  dashed  line  corresponds  to  the  magnetic 
insulation  threshold  for  a  PEOS .  The  solid 
line  corresponds  to  a  MITL  with  a  disruptiv 
structure. 
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low  voltage  region.  The  equations  governing  the  fits  are: 


1  )  z  /z 

1  p 


.51975  V 


.  15  + 


1  .51n(l  .  1  .95()7V)' 

1  +  1.9567V 


(164) 


2)  z , /z 

1  p 


.4675  V 


.  172 


.2 


1  + 


371n(l.  +  1.95G7V) 
1  +  1.95G7V 


.  {  165) 


where  equation  (164)  gives  the  magnetic  insulation  tlircshold 
for  simple  MITLs  and  equation  (165)  gives  the  threshold  for 
MITLs  with  a  PEGS.  The  voltage  difference  (V)  across  the 
anode/cathode  gap  is  in  megavolts.  These  formula  match  the 
simulation  data  to  within  a  percent.  They  arc  accurate  to 
within  2%  up  to  10  megavolts.  Extrapolating  out  to  10 
megavolts  is  made  possible  by  noting  the  similarities 
between  the  theoretical  predictions  of  Wang  [33]  and  the 
MASK  simulations  of  this  study.  It  should  be  a  simple 
matter  to  modify  the  above  relations  to  give  the  insulation 
threshold  out  to  over  a  100  roegavolts  using  the  same 
procedure.  The  magnetic  insulation  threshold  for 
unperturbed  MITLs,  as  given  by  equation  (164),  is  shown  in 
figure  38 . 

The  magnetic  insulation  threshold  for  a  system  with  a 
PEGS  (from  equation  (165))  is  shown  in  figure  39.  Also,  in 
figure  39.  is  a  plot  of  the  insulation  threshold  estimate 
found  and  used  at  the  Naval  Research  Laboratory.  Figure  40 
shows  the  corresponding  critical  currents  for  magnetic 
insulation  for  the  two  cases.  From  figures  39  and  40,  the 
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Figure  38.  Fit  to  Data  Corresponding  to  Figure  35. 

This  fit  is  for  the  cases  run  with  no  structures  and 
no  ion  emission. 
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Fit  to  Computational  Data  (PEGS,  Figure  37) 
NRL's  Previous  Estimate  (CRIT)  of  the  Magnet 
Threshold. 


Compared  to 
ic  Insulation 
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Figure  40-  Fit  to  Computational  Data  (ICNEW)  Compared  to  NRL's  Previous 
Estimate  (ICOLD)  of  the  Magnetic  Insulation  Threshold. 
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NRL  estimate  is  seen  to  be  very  close  to  that  obtained  in 
this  work  up  to  about  3  megavolts.  For  the  entire  voltage 
range,  the  current  required  for  magnetic  insulation  is 
greater  in  the  NRL  theory.  This  may  not  be  bad  in  practical 
applications  because  of  6  dependent  asymmetries  that  would 
degrade  the  magnetic  insulation.  To  insure  insulation, 
larger  values  of  the  critical  current  would  be  required.  In 
the  next  section,  equation  (165)  is  used  to  theoretically 
model  the  magnetic  insulation  and  opening  of  a  PEGS. 

5.4  Model  Gamble  II  Experiment 

In  section  5.3,  the  threshold  for  the  magnetic 
insulation  of  a  system  with  a  PEGS  is  given  (in  figure  39). 
To  verify  that  this  result  is  consistent  with  experiment, 
the  insulation  curve  of  figure  39  (equation  (165))  is  used 
in  conjunction  with  the  transmission  line  code  BERTHA  to 
model  a  Gamble  II  PEGS  experiment.  The  experimental  and 
theoretical  generator  and  load  currents  are  compared  along 
with  voltages  and  total  power  output.  Since  these 
quantities  are  closely  tied  to  the  magnetic  insulation 
phenomena,  a  close  match  between  the  theory  and  experiment 
will  justify  equation  (165). 

The  experimental  set  up  is  shown  in  figure  41  [31]. 

The  anode  is  at  a  radius  of  .05  meters.  The  cathode  radius 
is  .025  meters.  The  PEOS  plasma  is  produced  by  3  flashboard 
plasma  sources.  The  plasma  fills  a  .1  meter  long  region  of 
the  cathode/anode  gap.  The  diode  is  located  .05  meters 
downstream  of  the  PEOS  plasma.  The  plasma  source  fires  2 
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Figure  41.  Experimental  Setup  for  PEGS  Shot  3426. 


microseconds  before  the  Gamble  II  pulse  arrives  at  the 
switch.  This  is  to  allow  for  a  complete  plasma  prcfill 
[31].  The  data  for  an  opening  PEGS  (figure  15) 
corresponding  to  shot  3426  [31]  is  compared  to  the 
theoretical  BERTHA  predictions  using  the  new  insulation 
relation  (see  figure  39).  From  figure  42,  it  is  clear  that 
the  theory  is  doing  a  good  job  of  predicting  the 
experimental  results.  The  plasma  parameters  that  give  the 
properties  of  the  injected  plasma  correspond  to  those  of  the 
experiment.  However,  it  is  impossible  to  measure  the 
quantities  precisely.  As  a  result  there  is  latitude  for 
adjusting  them  to  improve  the  theoretical  prediction  given 
in  figure  42.  The  theoretical  predictions  here  are  similar 
to  those  obtained  by  Ottinger  et.  al.  in  previous  studies. 
This  can  be  understood  by  examining  figure  39  and  noting  how 
close  the  NRL  and  this  works  insulation  thresholds  are  in 
the  0  to  4  megavolt  regime  —  where  this  experiment  was 
car r i ed  out. 
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CHAPTER  6 
CONCLUSIONS 

In  Chapters  1  and  2  the  need  for  and  the  basic  theory 
governing  magnetic  insulation  in  coaxial  cylindrical  pulsed 
power  systems  are  introduced.  Some  theoretical  work 
previouly  done  in  cartesian  coordinates  is  extended  to 
cylindrical  coordinates.  In  comparing  the  resulting  theory 
to  MASK  simulations  it  is  found  that  the  cylindrictil  MITL 
electron  flow  is  laminar  in  nature,  but  has  a  non-Br i 1 1 ou  i  n 
f 1 ow  structure. 

This  motivates  the  derivation  of  a  general  relativistic 

laminar  MITL  flow  theory  in  Chapter  3  of  which  the  Brillouin 

flow  theory  is  a  special  case.  The  new  theory  is  used  to 

perform  an  independent  derivation  of  the  Brillouin  flow 

theory.  It  is  subsequently  used  to  generate  self-consistent 

MITL  flow  solutions  for  flows  with  constant  electron  density 

profiles.  Such  non-Br i 1 lou in  laminar  flows  are  found  to 

have  much  in  common  with  several  short  MITL  MASK 

simulations.  In  addition,  it  is  found  using  the  general 

theory  that  for  all  laminar  flows  (Brillouin  or  not)  the 

relativistic  factor  for  the  electrons  at  the  sheath  edge  is 

given  by:  t  =  I  /I  . 

mac 

In  Chapter  ■4,  the  electron  flow  downstream  of  abruptly 


varying  structures  (some  simulating  an  opening  PEGS)  is 
analyzed.  It  is  noted  that  such  structures  often  launch  a 
beam  of  magnetically  insulated  electrons  by  means  of  an 


23 


drift.  Then  the  regime  in  which  laminar  flows  downstream  of 
the  PEGS  occur  is  pointed  out.  A  set  of  equations,  obtained 
via  the  laminar  flow  theory  of  Chapter  3.  is  presented. 

These  equations  give  the  fields  and  properties  of  the 
launched  laminar  flow. 

In  Chapter  5,  the  magnetic  insulation  thresliolds  for 
coaxial  cylindrical  MITLs  are  given  for  MITLs  with  and 
without  disruptive  structures,  and  with  and  without  ion 
emission.  One  structure  of  particular  interest  simulates  a 
PEGS.  Variations  on  this  structure  lead  to  the  magnetic 
insulation  threshold  for  a  system  with  a  PEGS.  The  magnetic 
insulation  threshold  is  checked  for  the  PEGS  by  using  it  in 
a  transmission  line  code  to  predict  results  from  a  Gamble  II 
experiment.  The  theoretical  results  compare  well  with 
experimental  results.  This  corroborates  the  computationally 
determined  magnetic  insulation  threshold  for  the  PEGS. 

In  short,  this  work  presents  a  new  and  general  theory 
governing  laminar  flows  in  cylindrical  coordinates.  This 
theory  includes  the  Brillouin  flow  theory  as  a  special  case 
(corresponding  to  no  perturbations);  however,  the  general 
theory  can  generate  self  consistent  MITL  flow  solutions  for 
perturbed  systems.  This  work  also  sets  forth,  for  the  first 
time,  the  magnetic  insulation  thresholds  for  MITLs  with  and 
without  obstructions  and  with  and  without  ion  emission.  In 


addition,  the  magnetic  insulation  threshold  for  an  opening 
PEGS  is  given.  These  thresholds  are  given  simply  in  terms 
of  the  anode  and  cathode  currents,  geometry,  load  impedance. 
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vacuum  line  and  structure  impedance,  as  well  as  the  applied 
VO  1 tage . 

This  research  can  be  expanded  and  applied  to  various 
additional  ar'^as.  For  instance,  the  tecnniques  used  to 
obtain  the  general  theory  of  laminar  flows  in  cylindrical 
coordinates  could  be  used  to  extend  the  theory  to  other 
coordinate  systems.  One  could  also  use  the  general  theory 
to  obtain  the  magnetic  insulation  threshold  for  a  variety  of 
non-Br i 1 1 ou i n  flows.  The  general  theory  could  also  be  used 
to  generate  laminar  MITL  flows  for  pulsed  power  systems  with 
obstructions  and  impedance  mismatches.  In  addition,  the 
magnetic  insulation  threshold  curves  given  in  Chapter  5 
could  be  used  to  design  a  variety  of  pulsed  power  systems. 
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APPENDIX  A 

ALTERNATIVE  CYLINDRICAL  TRANSMISSION  LINE  THEORY  DERIVATION 

In  Appendix  A,  portions  of  the  cartesian  theory  of 
Mendel,  as  presented  in  reference  [10],  are  extended  to 
cylindrical  coordinates.  In  this  derivation,  the  function 
F{P^^.W)  needed  in  equation  (39)  arises  naturally.  This 
alternative  derivation  of  equations  (-43)  and  (44)  in  the 
text  is  included  in  order  to  add  additional  insights  to  the 
theory . 

This  derivation  differs  from  that  in  the  text  in  that 
no  Jacobian  of  transformation  is  explicitly  performed. 

Since  equations  (18)-(21),  and  (24)-(35)  still  apply  as 
shown  in  the  text  they  will  not  be  repeated  here.  What  is 
needed  are  the  electron  density  and  the  current 
distributions  defined  in  terms  of  the  canonical  momentum  and 
total  energy. 

The  electron  density  with  canonical  momentum  between 


P  and  P  +  dP  and  energy  between  W  and  W  +  dW  at  r  is 
c  z  c  z  c  z 

(from  conservation  of  charge): 


dP  dW  ■ 

cz  err 

where  |j(P  .W)|  is  the  one  way  current  density. 

C  Z 


(Al) 
The  factor 


of  2  comes  form  realizing  that  for  every  charge  entering  a 


unit  volume  from  the  top  there  is  one  also  entering  from  the 


bottom.  f  is  a  function  of  P^^  and  W  and  is  given  by  (see 
equation  (25) ) : 
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.2  2  .2  2  ,  .  ,2,n 

r  =  V  -z  =  V  -(c/tt)  (P^^+a) 

2 

Substituting  in  equation  (35)  and  eliminating  v  with 
equation  (3)  gives: 

r  =  .  (A2) 

Substituting  (A2)  into  (Al)  implies: 


2| j(P^^.W) |(4HW)dP^^dW 
ecr(((tHW)2-l-(P^^+a)2)^^2 


(A3) 


The  current  density  is  (A3)  times  z(r).  If  one  defines  2{r) 
in  terms  of  the  total  energy  and  canonical  momentum  (from 
equation  (25))  the  current  denstiy  is  given  by: 


dJ(Pc,.W) 


2(P^,^a)|j(P^^.w)|dP^^dW 

r((4KW)2_i-(P^^+a)2)^^2 


(A4) 


Now,  from  Poission’s  equation,  equation  (A3)  and  writing 
everything  d i mens i on  1 es s 1 y  gives: 


V 


2ep 

o 

me 


// 


1J(Pcz-^) 

r((ctH-W)2-i-(P^^+a)2)l/2 


(A5) 


From  Ampere’s  law,  the  restrictions  on  the  electron  flow 
given  in  section  2.1,  equation  (A4),  and  writing  the  result 
in  dimensionless  form  gives: 


V 


2 

a 

r 


2ep 

o 

me 


// 


r((4HW)2-l-(P^^+a)2)^^^ 


(A6) 
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Now,  from  (A5)  and  {A6)  the  electron  disrtibution  function 


is  de f i ned  to  be : 

2e/i 

F(P  .W)  =  - ■ 

^  C2  ^  me 


lj(P.,.W)| 


cz 


(A7) 


One  should  also  note  that  in  equations  {A5)  and  (A6)  that 
the  integrals  are  performed  in  the  canonical  momentum  and 
energy  phase  space.  The  variable  r  is  independent  of  this 
space  and  can  be  pulled  out  of  the  integrals.  So  from 
equations  {A5)-(A7)  and  the  definition  given  in  equation 
(40)  one  easily  obtains: 


-  i  . 

r  r  oa 

These  are  the  same  as  equations  (41)  and  (42)  which 
immediately  imply  equations  (43)  and  (44). 


(AS) 


(A9) 
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APPENDIX  B 

RELATION  BETWEEN  THE  DENSITY  AND  VELOCITY  PROFILES 
In  Appendix  B,  a  simple  yet  surprising  relation  is 
found  between  the  density  and  velocity  profiles  for  laminar 
flows.  This  relation  allows^one  to  obtain  a  complete 
self-consistent  laming  MITL  flow  solution  given  only  the 
density  profile,  or  magnetic  field  profile,  or  canonical 
momentum  and  total  energy  profile,  etc.  In  other  words,  a 
complete  flow  solution  is  possible  given  a  very  modest 
amount  of  information  about  the  flow. 

Deriving  the  relation  between  the  density  and  velocity 
profiles  begins  by  integrating  e^'a.TTon  (43 )  oT  Elie  TeXT 
from  the  cathode  to  some  r  less  than  or  equal  to  the  sheath 
radius  (r  ).  The  result  is: 


r(4)'2-a-2)  _  ^  gG  -  f  ( ({)’ ^-a  ’  ^ )  d  r  . 


(Bl) 


where  the  subcript  c  refers  to  the  cathode.  This  convention 
will  be  used  in  the  following.  The  subscripts  a  and  m  will 
refer  to  the  anode  and  sheath  edge  radius  respectively.  One 
should  also  note  that  G  =  0  for  laminar  flows  (equation  i 
( 40 ) )  because : 


c((4h-W)2-i_(p 

- - - ^ -  =  0  (B2) 


ADDENDUM: 


FJSRL-TR-89-0008,  Aug  89,  Magnetically  Insulated  Electron  Flows  in  Pulsed  Power 
Systems,  by  R.I.  Lawconnell  129 


Dimenslonalizing  (Bl)  implies 


E^Cr)  =  B^Cr^)  -  dr^ 

“  r  2 


(B3) 


r-r 


^  (c^B^Cr)  -  E^ir)) 


where  the  bar  signifies  an  average  over  r.  It  is  necessary 
to  write  B(r)  and  E(r)  in  the  following  form  using  Ampere’s 
and  Gauss’s  laws: 


B(r) 


2Trrd2 


It  I 
o  c 

2Trr 


ji  I 
'^o  c 


2~(VI,)f(r)g(r)  .  ^ 


E(r) 


9^(1') 

e  2Trrdz 
o 


(^a  - 

e  2Trrv  , 
o  d 


p ( r ) rdr 


(B4) 


(B5) 


where  r  and  r  are  the  cathode  and  sheath  radii 
c  m 

respectively.  p(r)  is  the  density  profile.  I  and  I  are 

3.  C 

the  anode  and  cathode  currents  respectively.  v^  is  the 
average  drift  velocity  and  equals: 


I  -I 
a  c 


J  "'p(r)27rrdi 


(B6) 
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q^Cr) 

dz 


enclosed  charge 
unit  length 


{B7) 


v^(r)  =  average  velocity  of  charge  with  radius  less 
or  equa 1  tor 


=  v^g(r) 


(B8) 


rr 


f{r)  = 


I  -I 
a  c 


p(r)2Trrdr  _  _ c 


[  p(r)rdr 
''  r 


'r 


(B9) 


m 


p{ r ) rd  r 


r 


p(r)v(r)2Trrdr 


g(r)  = 


a  c 


f  (r) 


(BIO) 


Now,  substituting  (B4)  and  {B5)  into  (B3)  yields 


_2,  .  2,_2  c„2,  ^ 

E  (r)  =  c  (B - -B^)  + 


1 


A  ^  f  ^2 

4ir  r  {  e  c  ) 

'  o  ' 


J 


r  lg-K2(I^-I^)I  Jg-(ce^27r)^E^{r)r^f^-^(I^-I^)(fg)^dr 


■  (Bll) 


For  simplicity  in  the  following  analysis,  this  is  evaluated 
between  the  cathode  (r  )  and  the  sheath  edge  radii  (r  ).  It 
becomes  after  algebra  (with  space  charge  limited  flow): 


2irr 


m 


2  2  1/2 

Si  c 


1+r  (kC  +1D  ) 
m '  m  m ' 


1  +  r  A 
m  m 


(B12) 
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where 


I  -I 

,  _  a  c 

''  =  r-Ti- 

a  c 


I  +I 
a  c 


-=  f"  ^  - 

r  r 
c 

D  =  r*" 


(B13) 


(BM) 


(BIS) 


(B16) 


(BIT) 


Now  for  laminar  flows,  from  equation  (50)  and  the  fact 
that  of  equation  (50)  equals  zero,  one  obtains  an 
additional  relation  for  the  electric  field  at  the  sheath 


edge  : 


^o  2  2 

E(r  )  =  15-2-  (I^  -  I^)  . 

m^  2irr  'a  c  •' 

m 


(BIS) 


Setting  equation  (BIS)  equal  to  equation  (B12)  and 
differentiating  with  respect  to  the  sheath  edge  radius  (t^) 
implies  after  algebra  that: 


1  g{r  ) 

f(rj  =  - -  .  (B19) 

1  -  kg2(r^) 

From  the  derivation  it  would  seem  this  relation  is  only 
valid  at  the  sheath  edge.  Nevertheless,  in  the  following  it 
will  be  shown  to  be  valid  not  only  at  the  sheath  edge,  but 
at  the  cathode  and  across  the  entire  electron  sheath.  Once 
this  relation  is  shown  to  be  true  across  the  electron  sheath 
in  will  be  used  to  relate  the  density  and  velocity  profiles. 
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First,  it  is  shown  that  {B19)  holds  at  the  sheath  edge. 

At  the  sheath  edge  (r  ).  f(r  )  =  1  =  g(r  )  from  equations 

'm  m  m 

(B9)  and  (BIO).  Substituting  g{r^)  =  1  into  equation  (819) 


yields: 


f(r  ) 
^  m 


1 


as  expec  ted . 

At  r  =  r^  ,  g(r^)  =  0.  Substituting  g(r^)  =  0  into 
equation  {B19)  yields  f{r^)  =  0.  At  r^ .  f{r^)  does  indeed 
equal  zero  as  can  be  seen  in  equation  (89). 

Now  it  is  necessary  to  show  that  for  laminar  flows 
(B19)  is  true  in  the  region  r^  <  r  <  r^.  For  laminar  flows: 

^  =  -V  X  I?  =  -v(r)B(r)  r  .  (820) 

Substituting  (84)  and  (85)  into  (820)  implies: 


f{r) 


I  v(r)v, 
c  '  d 


V  ,  1 


_  v(r)g(r)  d 

2 

c 


(821) 


where  is  given  in  this  case  by  equation  (52).  The 
velocity  v(r)  can  be  obtained  from  equation  (810)  assuming 
(B19)  is  true  across  the  flow: 


f{r) 


1  g(r) _ 

1  -  kg^(r) 


(822) 


Now  (822)  allows  one  to  write 
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g(r) 


2f(r)k 


[[ 


1 


2f{r)kJ 


The  derivative  of  (BIO)  gives: 


v(r) 


(g'  Cr)f  (r)  -«•  g(r)f  ■  (r))^^a 
2irrp  (  r  ) 


where 


f(r)  = 


g’  (r)  = 


The  re  f  ore , 


-  P(^)  ^ 


'I 


m 


p( r ) rdr 


I 

2kf^(r ) 


1  - 


1 


2kf(r)({ 


-J _ )2  ^  i) 

2f (r )k^ 


v(r)  = 


'J: 


2-?r|  p(r)rdr 
r 


[[ 


1  + 


{I-_I-)f-{r), 
1  1  a  c 

-1/2-1 

1  ,2  J 

. 

c 

Now  for  laminar  flows,  equation  (B21)  must  be 
satisfied.  If  it  is  satisfied  by  the  g{r)  and  v(r) 
from  assuming  (B22)  valid  across  the  flow,  then  the 
assumption  of  (B22)  is  justified.  So,  substituting 
(B25),  and  (52)  into  (B21)  yields: 


(B23) 


{B24) 


(B25) 


obtained 

{B23)  . 
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Given  an  anode  and  cathode  current,  and  f(r)  on  the 

right  hand  side  (RHS)  of  (B26).  the  calculated  f(r)  on  the 

left  hand  side  (LHS)  of  the  equation  must  equal  the  r(r) 

substituted  into  the  right  if  (B22)  is  to  be  justified. 

This  is  done  for  three  different  currents  and  two  different 

f(r)s  in  Table  Bl.  Equation  (B26)  not  only  holds  for  the 

cases  given  in  Table  Bl,  it  holds  in  general.  This 

justifies  equation  (B22)  across  the  flow.  Therefore. 

f{r)  =  I -  .  r  <  r  <  r  .  (B27) 

1  -  kg  (r) 

Equation  (B27)  is  used  in  the  text  to  determine  a  velocity 
profile  given  a  density  profile. 


f {r)RHS 

I 

a 

I 

c 

g(r) 

Num . 

Den  . 

f ( r )LHS 

.  5 

I 

.91 

.5202592  . 

4859543 

. 97190S7 

.  5 

a 

a 

.5 

I 

.51 

.6457513 

3779645 

.755929 

.5 

a 

a 

.  25 

I 

.31 

.4756728 

1956984 

.7827937 

.  25 

a _ 

a 

Tab  1 e  B 1  . 

Verification  of  equation 

(B26)  . 

This  is 

done  for 

three  different  anode/cathode  current  ratios  and  two 
different  f{r)  values. 


APPENDIX  C 


MAGNETIC  INSULATION  THRESHOLD  DATA 
Tables  Cl  thru  C9  comprise  a  record  of  pertinent 
magnetic  insulation  data  for  100  MASK  simulations.  The 
simulations  were  run  for  coaxial  cylindrical  MITL 
geometries,  with  and  without  structures,  and  with  and 
without  ion  emission  off  the  structures.  The  ion  emitting 
structures  were  used  to  duplicate  PEOSs  in  different  stages 
of  opening.  Diagrams  of  the  two  fundamental  computational 
meshes  used  in  this  study  are  given  in  figures  3a  and  3b. 

The  mesh  parameters  in  3a  and  3b  are  varied  to  obtain  an 
understanding  of  the  relation  between  the  applied  voltage, 
geometry  and  magnetic  insulation.  Tables  C1-C3  correspond 
to  figure  3a  and  Tables  C*l-C9  correspond  to  figure  3b.  The 
data  from  these  tables  is  used  to  generate  the  magnetic 
insulation  threshold  curves  for  figures  35  and  37  in  Chapter 
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S i mu  1  a  t i on 

V 

d 

I 

I 

Z 

a 

g 

a 

c 

0 

M407-M412 

2x10^ 

.01 

-1 .56x10® 

-1.5x10® 

.31 

M414 

” 

.015 

-9.6x10^ 

o 

X 

00 

1 

.  o 

M416 

.  032 

-7.92x10"^ 

-4.72x10^ 

.  61 

M417.M419 

»« 

.  023 

-7.872x10^ 

-4.27x10^^ 

.61 

M418.M420 

«* 

.018 

-8,46x10^ 

-6.5x10^ 

.  568 

M421 

f « 

.016 

-9.3x10® 

— 

.  52 

M423 

«« 

.01  1 

-1 . 324x10® 

— 

.  363 

M424 

•• 

.02 

-7.97x10^ 

— 

.  603 

M426 

.5x10^ 

.015 

-2.753x10^ 

-1 . 79x10^^ 

.  437 

M427 

•  t 

.011 

-2.784x10"^ 

-1 .87x10"^ 

.  431 

M435 

o 

X 

.016 

-4.65x10^ 

-2.672x10^ 

.517 

M436 

«« 

.  018 

-4.613x10^ 

-2.463x10'^ 

.521 

M437 

It 

.014 

-4.846x10"^ 

-2.707x10"^ 

.  496 

M438 

It 

.012 

-5. 169x10^ 

— 

.  465 

M440 

3x10® 

.  022 

-1.118x10® 

-4.446x10"^ 

.  645 

M441 

•• 

.024 

-1 .095x10® 

-4  .  ISSxlo'^ 

.  659 

M442 

•• 

.020 

-1 . 154x10® 

-7.901x10® 

.  625 

M443 

- 

.  028 

-1 .093x10® 

-4  .  163x  lO'^ 

.  660 

Table  Cl  Simulation  data  for  the  determination  of  the 

magnetic  insulation  threshold  (r  =  .05.  r  =  .025).  This 

St  c 

study  was  done  without  ion  emission  off  the  anode. 
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Simulation 

V 

d 

I 

I 

z 

a 

g 

a 

c 

o 

M450 

2x10® 

.056 

-2.509xl0'^ 

-1 . 203x10^ 

.  577 

M451 

«« 

.030 

-2. 481x10^^ 

-1 . 136x10^ 

.  5835 

M452 

•• 

.012 

-2.623x10"^ 

-2.  Ixio'^ 

.  552 

M453 

.020 

-2.44x10^ 

-1.11x10^ 

.  593 

M456 

«t 

.018 

-2.45x10"^ 

.  591 

M457 

1x10® 

.022 

-1 .44x10^ 

-1 . 44xl0'^ 

.  501 

M458 

•  t 

.006 

-1 .77x10^^ 

-1 . 77x10^^ 

.  408 

M459 

•• 

.016 

-1 .43x10^^ 

-1 . 43x10^ 

.  505 

Table  C2 
magne  tic 
study  is 


Simulation  data  for  the 
insulation  threshold  (r 

a 

without  ion  emission  off 


determination  of 
=  .05.  r^  =.005) . 
the  anode . 


the 
Th  i  s 
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Simulation 

V 

d 

I 

I 

z 

a 

g 

a 

c 

o 

MNOl 

2x10® 

.01 

-1.779x10® 

-1 .365x10® 

.  557 

MN04 

•> 

.018 

-1.646x10® 

-7 . 37x10^^ 

.602 

MN05 

- 

.014  - 

-1 .642x10® 

-7 . 305x10^ 

.  6033 

MI30 

2x10® 

..014  - 

-1 .804x10® 

— 

.  549 

MI31 

tt 

.024  - 

-1 .737x10® 

— 

.  5703 

MI33 

t* 

.018  - 

-1 .798x10® 

— 

.  551 

MI35 

3x10® 

.024  - 

-2.518x10® 

— 

.  59 

MI36 

«* 

.018  - 

-2.491x10® 

— 

.  5965 

MI38 

.014  - 

-2.674x10® 

— 

.55 

Table  C3 .  Simulation  data  for  the  determination  of  the 
magnetic  insulation  threshold  (r  =  .035,  r  =  .025).  The 

3i  C 

"MN"  calculations  were  done  without  ion  emission  and  the 
"MI"  simulations  were  done  with  ion  emission  off  the  anode. 
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Simulation 

V 

a 

d 

g 

I 

a 

h 

z 

P 

Flow 

PN05 

2x10® 

.0075  - 

-1 .834x10® 

.54 

vortex 

PN04 

- 

.008 

-1 .898x10® 

.52 

vortex 

PN06 

t« 

.009 

-1 .509x10® 

.656 

1  am i nar 

PN07 

•  1 

.012 

-1 . 473x10® 

.  672 

1  am i nar 

PN08 

•• 

.015 

-1 .469x10® 

.674 

1  ami nar 

PN09 

•  • 

.0105  - 

-1 . 488x10® 

.666 

1  am  i  nar 

PNIO 

3x10® 

.0105  - 

-2.025x10® 

.734 

1  ami nar 

PNl  1 

- 

.012 

-2.023x10® 

.735 

1  am  i  nar 

PN12 

•• 

.0075  - 

-2.63x10® 

.565 

vortex 

PNl  3 

If 

.021 

-2.093x10® 

.71 

1  am i nar 

PN14 

4x10® 

.012 

-2.569x10® 

.7713 

laminar 

PN15 

f  f 

.0195  - 

-2.633x10® 

.7525 

laminar 

PN16 

11 

.009 

■2.858x10® 

.6933 

vortex 

Table  C4 .  Simulation  data  for  the  determination  of  the 

magnetic  insulation  threshold  for  systems  with  a  non-ion 

emitting  structurefr  =  .025.  r  =  .05,  r  =  .035,  D  =  . 

c  a  p 

D  =  .06.  Az  =  .03) . 

P 


192, 
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Simulation 

V 

a 

d 

g 

I 

a 

z 

P 

flow 

PN12 

4x10® 

.0075 

-3.615x10® 

.548 

vortex 

PN18 

'• 

.012 

-2.573x10® 

.7701 

1  am i nar 

PN19 

.0195 

-2.639x10® 

.  7501 

•• 

PN21 

3x10® 

.009 

-2.074x10® 

.717 

1  am i nar 

PN22 

ft 

.0105 

-2.023x10® 

.  735 

PN23 

ft 

.0135 

-2.06x10® 

.7214 

tf 

PN24 

2x10® 

.0105 

-1 .476x10® 

.671 

1  am i na  r 

PN25 

•• 

.009 

-1 . 491x10® 

.665 

•• 

PN26 

ft 

.0075 

-1 .829x10® 

.542 

V o r  t  ex 

PN27 

6x10® 

.0106 

-3.56x10® 

.835 

laminar 

PN28 

f  f 

.015 

-3.623x10® 

.82 

•  1 

PN29 

tf 

.0225 

-3.774x10® 

.788 

•• 

Table  C5 .  Simulation  data  for  the  determination  of  the 

magnetic  insulation  threshold  for  systems  with  a  non-ion 

emitting  structure  (r  =  .025,  r  =  .05,  r  =  .035,  D  = 

c  a  p  p 


.105,  D 


.192,  Az 


.  03)  . 


Ml 


Simulation 

V 

a 

d 

g 

I 

a 

z 

P 

Flow 

pool 

2x10^ 

.01 

-1 .533x10® 

.  463 

vortex 

M482 

«« 

.016 

-1 .067x10® 

.665 

1  ami na  r 

M483 

•  1 

.02 

-1 .071x10® 

.662 

Where 

r  = 

.  04, 

D  =  .07. 

D  =  .256.  Az  =  .05 

P 

P 

PN31 

4x10® 

.018 

-1 . 83x10® 

.  7751 

laminar 

PN32 

.012 

-2.492x10® 

.  569 

vortex 

PN33 

- 

.024 

-1 .846x10® 

.  768 

1  am i na  r 

Where 

r  = 

.04. 

D  =  .06. 

D  =  . 192.  Az  =  .03 

P 

P 

PN35 

4x10® 

.018 

-1 .753x10® 

.647 

vortex 

PN36 

.012 

-2.525x10® 

.  449 

vortex 

PN37 

*• 

.024 

-1.469x10® 

.772 

laminar 

PN38 

.030 

-1.48x10® 

.766 

M 

Where 

r  = 

.045. 

D  =06. 

D  =  .192,  Az  =  .03 

2 _ E 


Table  C6 .  Simulation  data  for  the  determination  of  the 
magnetic  insulation  threshold  for  systems  with  a  non-ion 
emitting  structure  (r  =  .025,  r  =  .05). 


M2 


S i mu  1  a  t i on 

V 

a 

d  I 

g  a 

z 

P 

Flow 

PI  03 

2x10^ 

0075  -1.868x10^ 

.  53 

1  am i na  r 

PI04 

009  -1.844x10^ 

.5372 

•• 

PI05 

It 

0105  -1.84x10^ 

.538 

•  1 

PI06 

•• 

0135  -1.72x10^ 

.  575 

PI  08 

3x10^ 

0075  -2.664x10^ 

.  558 

1  am  i  na  r 

PI09 

>• 

009  -2.53x10^ 

.  587 

PI  10 

11 

0105  -2.503x10^ 

.  594 

•• 

PI13 

4x10^ 

0075  -3.492x10^ 

.  567 

V  o  r  t  cx 

PI14 

-• 

009  -3.222x10^ 

.615 

lam i na  r 

PI15 

•• 

0105  -3.165x10^ 

.  626 

PI  16 

11 

012  -3.148x10^ 

.629 

Table 

C7. 

Magne  tic 

insulation  threshold  data 

for  an  ion 

emitting  PEGS  like 

structure  (r  =  . 

'  c 

025.  r  = 
a 

.05,  T  = 

P 

.035. 

D  = 
P 

.06.  D  = 

. 1 92 .  and  Az  =  . 

03)  . 

1  43 


S i mu  1  a  t i on 

V 

a 

d 

.  g 

I 

a 

z 

P 

Flow 

PI  17 

3x10® 

.0075 

-2.567x10® 

.58 

1  am i na  r 

PI  18 

•• 

.  009 

-2.498x10® 

.  595 

PI  19 

•• 

.0105 

-2.491x10® 

.  5966 

PI21 

4x10® 

.009 

-3. 165x10® 

.626 

1  am i na  r 

PI22 

" 

.012 

-3. 151x10® 

.629 

•• 

PI23 

«« 

.0135 

-3. 13x10® 

.633 

PI26 

6x10® 

.012 

-4.376x10® 

.679 

laminar 

PI27 

*« 

.009 

-4.454x10® 

.667 

•• 

Table  C8 .  Magnetic  insulation  threshold 

emitting  PEGS  like  structure  {r^  =  .025, 

.035,  D  =  .09.  D  =  .192,  and  Az  =  .03). 
P 


data  for  an 

r  =  . 05  ,  r 
a  p 


i  on 
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S i mu  1 a  t i on 

V 

d 

I 

z 

flow 

a 

e 

a 

P 

PI29 

4x10® 

.021 

-2.247x10® 

.631 

laminar 

PI30 

ft 

.015 

-2.244x10® 

.632 

•• 

PI31 

•• 

.012 

-2.244x10® 

.632 

•• 

PI34 

6x10® 

.012 

-3. 145x10® 

.677 

1  am i na  r 

PI36 

ft 

.  027 

-3. 128x10® 

.68 

Table  C9 .  Magnetic  insulation  threshold  data  for  an 
emitting  PEGS  like  structure  (r  =  .025,  r  =  .05,  r 
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